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FORKWORD 


Initial  auti-ority  to  ;cnduet  5.  model  investigation  for  the  Hilo 
Hari-or  Tsunami  ana  navigation  Project  was  granted  by  the  Office,  Chief  of  . 
Engineers,  on  r~  Bept ember  19&1 .  The  U.  S.  Army  Engineer  Waterways  Experi¬ 
ment  Sta.tion  (WES)  was  authorised  to  resign  a  wave  machine  for  the  Hilo 
Harbor  tsunami  model  in  a  teletype  received  1 6  August  19^3  from  the 
U.  8,  army  Engineer  District,  Honolulu. 

This  analytical  investigation  was  conducted  by  Dr.  G.  H.  Keulegan , 
consultant,  ViEo  Hydraulics  Division,  simultaneously  with  tests  to  investi¬ 
gate  design -of;  a  bore  generator  for  the  Kilo  Harbor  model,  during  the 
period  November  19b3-May  19o4,  in  the  Water  Waves  Branch,  Hydraulics  Divi¬ 
sion,  WES,  under  the  direction  of  Mr.  E.  P.  Fort  son,  Jr.,  Chief  of  the 
Hydraulics  Division,  and  Mr.  K.  Y.  Hudson,  Chief  of  the  Water  Waves  Branch. 
The  Fortran  program  for  the  numerical  solution  of  equations  93?  9^,  and  95 
was  worked  out  by  Mr.  Michael  Dorl  of  the  Hydraulic  Analysis  Branch.  This 
report  was  prejared  by  Dr.  Keulegan. 

Directors  of  the  WES  luring  the  conduct  of  this  investigation  and 
preparation  of  this  report  were  Col.  Alex  G.  Sutton,  Jr.,  CE,  and 
Col.  John  K.  Oswalt,  Jr.,  C-S.  Technical  Director  was  Mr.  J.  B.  Tiffany. 
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t>  ASSTAACT 

After  a  brief  discussion  of  the  genesis  of  oscillatory  waves  by  the 
pneumatic  method,  the  required  pressure  condition  for  the  genesis  of  long 
surges  with  constant  wave  heights  is^.  inquired-into .  Reasoning  frcta  the 
energy  point  of  view,  it  is  seen  that  the  sum  of  the  pressure  head  in  the 
pneumatic  chamber  and  the  water  elevation  therein  should  remain  constant 
during. the  issuance  of  the  /ave  from  the  generator. 

Experience  indicates  that  this  condition  would  be  realized  for  the 
operation  where,  subsequent  to  the  raising  of  water  into  the  chamber,  one 
allows  the  outside  air  to  enter  the  chamber  through  an  aperture.  Because 
of  the  inertia  of  water,  marked  pressure  oscillations  occur  in  the  chamber 
at  the  instant  cf  aperture  opening.  These  oscillations  last  for  a  short 
time,  and  scon  the  fall  of  water  surface  and  the  increase  of  pressure, 
Nboth,  are  uniform.  These  aspects  of  the  generator  behavior  are  amenable 
tp  analysis. 

The  design  of  a  pneumatic  generator  for  long  surges  may  be  effected 
solely  on  the  basis  of  latter  manifestations  of  the  generator  action.  ‘  Pro¬ 
cedures  for  the  design  of  a  pneumatic  generator  for  long  waves  are  outlined 
and  are  applied  for  the  generator  used  in  the  Hilo  Harbor  model  used  to 
study  the  problems  of  protection  against  tsunamis.  . 
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GLOSSARY 


Pneumatic  generator  and  channel 
D  Depth  of  nozzle  tliroat 
h  Wave  height,  surge  or  oscillatory  wave 
h  Maximum  wave  height  in  oscillatory  wave 
H  Maximum  height  the  water  is  raised  by  suction  in  chamber 
Hg  Base  height  of  the  prismatic  portion  of  channel 
HQ  Depth  of  water  in  channel  downstream  of  wave  front 

Height  of  generator  ceiling  measured  from  channel  bottom 
k  Wave  number;  k  =  2if/X 
Z  Length  of  pneumatic  chamber 
Z^  Length  of  nozzle 

Q  Volume  or  area  of  oscillatory  wave  in  channel  of  unit  width 
positive  or  negative  portion 

t  Time 

tg  Characteristic  time;  See  equation  8l 
T  Period 

V  Value  of  air  gap  during  wave  generation;  V  =  i(5  +  a)  •  1 
VQ  Initial  value  of  air  gap;  VQ  =  iA  •  1 
6  Fall  of  water  surface  in  chamber  during  wave  generation 
A  Initial  air  gap  in  generator;  =  H  +  A 


t-j  •  j  ^  ^  .  _  A  •>  An 

i  GHuu  uuiuut?x  ^  u  —  c:1'  j  -L 


Velocities 


2  2  2 
q  Absolute  velocity;  q  =  u  +  v 

u  Horizontal  component  of  velocity  or  particle  velocity  in  wave 
motion 

u  Maximum  particle  velocity  in  oscillatory  wave 
m 

Ug  Particle  velocity  in  s’irges  from  generator 
v  /ordeal  component  of  velocity 
v  Velocity  in  the  direction  of  the  normal  to  a  surface 
v.  Fall  velocity  of  water  surface  in  chamber 
-o  Wave  velocity  in  surges 


Pressures 


A  Semiamplitude  of  the  pressure  variations 

k, ,k  Constants  to  describe  the  fall  of  pressure  in  the  pneumatic 
*  2  chamber  for  surges  or  long  waves;  Ap  =  (l  -  k^  -  2k^t )ApQ 

pQ  Atmospheric  pressure 

P  Pressure  in  chamber  measured  relative  to  atmospheric  pressure 

PQ  Pressure  in  water  at  nozzle  mouth  measured  relative  to  atmospheric 
pressure 

Ap  air  suction  pressure  in  chamber  during  wave  generation 

Ap.  Air  suction  pressure  in  chamber  during  the  start  of  waves 

Ap^  A  positive  quantity;  represents  the  additional  pressure  difference 
to  overcome  the  friction  resistance  through  the  air  nozzle 

A',  suction  pressure  to  raise  water  to  elevation  H  ; 

U  Ap0  =  Pg(H  -  Kq) 

Energies 


E  Total  energy;  E  =  T  +  V 


viii 


Ej,  Internal  kinetic  energy  of  air  in  the  vessel 

«% 

E,  Rate  of  dissipation  of  kinetic  energy  from  the  surface  of  contact 
g  Constant  of  gravity 

Ie  An  expression  equal  to  - J*  “  l^^p- 

M  A  numerical  factor  in  the  e  tpressica  relating  the  kinetic  energy 
of  liquid  in  the  chamber  and  the  nozzle  to  h  and  HL  ;  see 
equation  33 

M  A  numerical  factor  in  the  expression  relating  the  kinetic  energy 
of  liquid  in  the  chamber  and  the  nozzle  xo  v1  and  i  ;  see 
equation  86 

A  numerical  factor  in  the  expression  relating  the  kinetic  energy 
of  liquid  in  the  nozzle  to  Ug  and  Ky  ;  see  equation  2.6 

T  Kinetic  energy  of  moving  liquids 

Tjj  Kinetic  energy  in  the  lower  part  of  the  pneumatic  chamber 

T2  Total  energy,  T2  =  T21  +  Tgg 

T21  Kinetic  energy  in  J  pneumatic  chamber  and  the  nozzle  below 
Tgg  Kinetic  energy  in  the  wave  in  the  channel 
V  Potential  energy  of  liquids 
W  Work  done  by  pressures 

£E  Difference  in  energies' of  a  portion  of  liquid  for  the  instants 
t1  and  t2 

a  Gravitational  potential:  n  =  fey 
Plow  of  air  into  receiving  vessels 

2 

a  Cross-sectional  area  of  orifice;  a  =  «d  /k 

Cq  Velocity  of  sound  in  the  air  outside  vessel 

Cq,C  Initial  and  later  discharge  coefficients,  respectively,  real  or 
J  apparent,  of  orifice 

Cp  Specific  heat  at  constant  pressure 

C  Specific  heat  at  constant  volume;  C  =  C  +  R 
v  p  v 


ix 


Diameter  of  air  orifice 


K 


Kiunerieal  defined  by 


m  Mass  of  gas  in  receiving  chamber 

Momenta *  of  liquid  entering  the  orifice 
n  A  constant 

Np  numerical  defined  by  =  7Pq/^Pq 

Pq  Pressure  of  outside  air 

p^  Pressure  at  vena  contracta  of  the  orifice 

Pressure  inside  receiving  vessel 

q1  Velocity  of  air  at  the  vena  contracta 

R  Gas  constant;  pv  =  RB 

v  Specific  volume;  v  =  l/p 

V.  Volume  occupied  by  air  jet 

VQ  Initial  volume  of  receiving  chamoer 

X.  Length  of  air  jet 
d 

6  A  positive  constant  _ 

7  Ratio  of  specific  heats;  7  =  C ^/C.^ 

A  Initial  ai*-  gap  between  the  water  surface  and  the  chamber  ceiling 

e  Internal  energy  of  gas  per  unit  mass;  or  a  numerical  which  depends 
on  the  circumstances  of  the  filling  and  which  can  be  a  variable 
changing  with  time 

6  Absolute  temperature 

t?  Temperature  in  degrees  Rankine  corresponding  to  32  F 
u 

p  Density 

pQ  Density  of  air  outside 

p.,  Density  of  air  ac  vena  contracta 


x 


.  ?.T^.’-3*.  tv/  — 


Density  of  air  inside 
Cross-sectional  area  of  vena  contracta 
GeoEetrical  ^aantities 

£,m  Directional  cosines  of  normal  drawn  inward  of  S 
s  Length  of  curve,  closed  or  open 
S  Surface  area 
x  Longitudinal  coordinate 
y,z  Vertical  coordinates 
Dimensionless  parameters 

”  1+*fc 


£  Hozzle  coefficient  of  resistance;  see  equation  143 

i  hA) 

9  6/A 

X  Hozzle  factor  of  resistance;  see  equation  52 
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NUMMARY 


After  a  brief  discussion  of  the  genesis  cf  oscillatory  waves  by  the 
pneumat!  *  method,  the  required  pressure  condition  for  the  genesis  of  long 
surges  with  constant  wave  heights  is  inquired  into.  Reasoning  from  the 
energy  point  of  view,  it  is  seen  that  the  sum  of  the  pressure  head  in  the 
pneumatic  chamber  and  the  water  elevation  therein  should  remain  constant 
during  the  issuance  of  the  wave  from  the  generator. 

Experience  indicates  that  this  condition  would  be  realized  for  the 
operation  where,  subsequent  to  the  raising  of  water  into  tne  chamber,  one 
allows  the  outside  air  to  enter  the  chamber  through  an  aperture.  Because 
of  the  inertia  of  water,  marked  pressure  oscillations  occur  in  the  chamber 
at  the  instant  of  aperture  opening.  These  oscillations  last  for  a  short 
time,  and  soon  the  fall  of  water  surface  and  the  increase  of  pressure, 
both,  are  uniform.  These  aspects  of  the  generator  behavior  are  amenable 
to  analysis. 

The  design  of  a  pneumatic  generator  f<  •  long  surges  may  be  effected 
-solely  on  the  basis  of  lattei  manifestations  of  the  generator  action.  Pro¬ 
cedures  for  the  design  of  a  pneumatic  generator  for  long  waves  are  outlined 
and  are  applied  for  the  generator  used  in  the  Hilo  Harbor  model  used  to 
study  the  problems  of  protection  against  tsunamis. 


APPROXIMATE  THEORIES  OP  PNEUMATIC  WAVE  GENERATORS 


/draulic  Laboratory, 


PART  I:  INTRODUCTION 


1.  The  first  phase  of  the  Hilo  Bay  model  investigation  will  comprise 
examination  of  the  effectiveness  of  various  protective  structures  against 
tsunamis  using  a  singular  positive  wave  of  translation  of  great  length  with 
large  and  uniform  heights  entering  the  bay.  Such  a  wave  can  be  created  in 
a  model  by  suddenly  releasing  the  impounded  waters  o*  a  reservoir,  by  uni¬ 
formly  displacing  a  bulkhead  through  a  distance,  or  by  using  a  pneumatic 
wave  generator  of  considerable  capacity.  The  elementary  aspects  of  the 
theory  relating  to  the  latter  will  be  discussed  herein. 

2.  The  original  idea  of  a  pneumatic  wave  generator  goes  back  to  the 

late  Professor  R.  T.  Knapp,  California  Institute  of  Technology,  who  used 

such  a  machine,  in  a  model  to  study  the  long  wave  conditions  of  Apra  Harbor, 

Guam.  Professor  Knapp  was  well  known  for  his  ingenious  and  original  ideas 

in  mechanical  devices,  and  the  air-operated  machine  unquestionably  is  a 

1* 

firm  testimonial  to  his  inventive  genius.  In  an  unpublished  paper  Knapp 
briefly  touched  upon  the  theory  of  the  machine  that  was  used  for  the  Apra 
Harbor  model.  In  kindly  transmitting  a  copy  of  the  paper  to  this  writer, 
Mr.  J .  M.  Caldwell  of  the  U.  S.  Coastal  Engineering  Research  Center  added 
remarks  to  the  effect  that  there  is  difficulty  in  following  the  arguments 
of  the  paper.  The  brevity  of  the  treatment  and  possible  typographical 
errors  could  be  the  reasons  for  this  difficulty1.  The  explanations  of  Knapp 
purport  to  deal  with  the  machines  for  the  purpose  of  generating  oscillatory 
translation  waves.  In  the  following,  before  entering  into  the  discussion 
of  air  machines  to  generate  singular  waves,  it  would  be  .instructive  to  con¬ 
sider  first  the  matter  of  the  generator  for  long  oscillatory  waves.  This 
is  done  for  the  purpose  of  developing  a  feeling  about  the  all-round  aspects 
of  the  problem. 

*  Raised  numerals  refer  to  similarly  numbered  items  in  the  Literature 
Cited  at  end  of  text. 
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'.-/hen  a  pneumatic  generator  is  provided  with  an  exit  nozzle,  and 
this  seems  to  !>'■  necessary  for  the  genesis  of  oscillatory  shallow-water 
waves  of  long  length,  the  flow  of  water  from  the  pneumatic  chamber  into  the 
channel  outside  is  one  of  deflection.  Owing  to  the  latter  condition,  it 
would  appear  that  the  mechanics  of  flow  are  more  readily  treated  by  resort¬ 
ing  to  energy  considerations  rather  than  to  the  momentum  principle.  The 
energy  relation  can  be  expressed  in -two  ways,  one  representing  the 
Lagrangian  point  of  view  and  the  other  the  Eulerian  point  of  view.  In  the 
review  of  Knapp's  problem,  both  forms  of  the  energy  relation  were  used. 

This  was  successful,  and  accordingly  the  same  approach  was  incorporated  in 
the  problem  of  a  pneumatic  wave  generator  to  produce  a  singular  impulse  of 
wave  of  great  length  and  of  constant  height  moving  over  water  of  constant 
depth  and  width. 

4.  The  present  report  is  not  meant  to  be  a  systematic  treatment  of 
pneumatic  generators.  It  merely  presents  the  original  thinking  that  guided 
the  conduct  of  experiments  on  a  few  shapes  of  generators.  The  order  of  the 
material  shown  here  does  represent  the  progress  of  analysis  running  parallel 
v/ith  the  tests.  The  aim  was  to  develop  certain  basic  relations  that  by 
themselves  were  sufficient  for  design  of  the  generator  which  could  be 

used  in  the  Hilo  Bay  model.  These  relatic  s  are  shown  in  the  section 
dealing  with  design. 

5.  For  the  interpretation  of  the  various  theoretical  results  in  the 

paper,  reference  is  made  to  some  of  the  experimental  results  with  various 

shapes  of  pneumatic  generators.  Test  procedures  leading  to  these  data  are 

2 

fully  explained  by  Mr.  C.  C.  Shen  in  a  separate  report. 
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FART  II:  THEORY  OF  FNFUMATIG  GENERATORS  FOR  OSCILLATORY  WAVES 


Kinematic  Relations 


-  6.  Definitions  of  symbols  for  -wave  height  and  particle  velocity 
(translation  waves)  are  as  follows: 

h  Wave  elevation  measured  from  undisturbed  surface 
Semiamplitude  of  wave  height,  maximum  value 
Hq  Depth  of  undisturbed  water  in  channel 
T  Period 

u  Uniform  particle  velocity  In  a  section 
u  Semiamplitude  of  particle  velocity,  maximum  value 
x  Distance  measured  from  moutn  of  nozzle 
X  Wave  length 

The  mouth  of  the  nozzle  is  that  opening  where  the  mean  depth  of  water  dur¬ 
ing  a  complete  oscillation  equals  HQ  .  It  is  assumeu  that 


and 


h  =  h  sin  (at  -  kx);  a  =  2«/r 

Q 

u  =  u^  sin  (crt  -  kx);  k  =  2«/R 


The  relation  between  h  and  u  is  determined  from  the  condition  c  f 

m  a 

continuity 

ah  cos  (at  -  kx)  -  u  HJc  cos  (at  -  kx)  =  0 
m  m  O  ’ 

„  h 
a  m 


or 


u 

m 


•V  * 

A  _m 

T  IL 


(i) 
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7.  Lot  Q  bo  the  area  of  the  elevated  portion  of  the  wave  or  the 
depressed  portion  in  channel  of  unit  width.  See  fig.  la.  For  t  =  0 


h  =  h  s:*n  kx 
m 


Q  =  %  f  sin  kx  d> 

J>./2 


2h 


Q  = 


m 


h  a 
m 


(2) 


k  rt 

8.  Under  the  action  of  sinusoidally  varying  pressure  P  ,  the  sur¬ 
face  of  water  in  the  pneumatic  chamber  during  a  complete  cycle  is  forced 
from  level  A  to  level  B  and  then  is  pulled  back  to  the  initial  level  A 
The  average  level  is  in  the  same  plane  as  the  undisturbed  water  surface  in 
the  outside  channel.  Let  the  displacement  5  of  the  water  surface  in  the 
chamber  be  measured  from  the  undisturbed  level.  By  storage  conditions, 
where  i  is  the  chamber  length. 


i  =  -uH 
dt  0 


where  u  is  now  the  particle  velocity  at  the  mouth  of  the  nozzle,  that  is 

at  :<_=  0  .  As  u  is  sinusoidal,  assume  that  6  is  also  sinusoidal  and, 

as  5  varies  from  5  to  -6  ,  then 

m  m 


5=5  cos  at 
m 

At  the  nozzle  mouth,  that  is,  at  section  x  ~  0  ,  the  particle  velocities 
are 


u  =  u  sin  oa 
m 


This  and  the  storage  equation  yield 


2*5 

a 


u  HnT 
m  0 


h  X 


k 


or 


(3) 


2lb  =  Q. 
m 

Accordingly,  the  pneumatic  chamber  will  require  a  storage  of  only  Q  . 

9.  To  produce  waves  of  specified  wave  heights  2hffl  ,  the  pressure  P 
of  the  pneumatic  chamber  needs  to  be  properly  controlled.  To  determine  the 
relation  of  P  to  h  ,  resort  may  be  made  to  the  equation  of  energy  to  be 

-IQ 

discussed  below. 


Energy  Equations 

10.  For  the  case  of  two-dimensional  flow,  the  equations  of  motion  are 
given  as 


du  1  dp  do 

dt  ~  p  dx  ox 

dv  =  _1  d£  _  dO 

dt  ”p  dy  <5y 


(*> 


where  d/dt  denotes  differentiation  along  a  particle  path,  that  is 


d  d  d  d 

at  ”  St  +  "  3k  +  v  3? 

where  u  and  v  are  the  velocity  components  in  the  coordinate  axis  x  and 

y  .  The  physical  meaning  of  ft  is  that  it  denotes  the  potential  energy, 
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per  unit  mass,  at  the  point  x,y  in  respect  to  gravitational  forces.  As 
both  p  ,  the  pressure,  and  ft  ,  the  gravitational  potential,  enter  into 
the  equation  as  rates,  the  resulting  motions  are  not  altered  if  the  value 
of  p  is  reduced  by  one  constant  and  ft  by  another  c/ustant.  Tills  allows 
us,  when  dealing  with  gravitational  wave  motion ;  to  put  atmospheric  pres¬ 
sure  equal  to  zero.  The  mathematical  expression  for  ft  is  gy  if  y  is 
drawn  vertical.  Now  y  may  be  measured  from  any  horizontal  plane.  For 
the  problems  ct  riand,  y  will  be  measured  from  the  channel  bottom. 

11.  The  derivation  of  the  energy  equation  on  the  basis  of  flow 

3 

equations  is  given  by  Lamb.  For  incompressible  flow  in  two 


5 


d:mensio as,  the  equation  reduces  to 


d 

at 


(T  +  V)  =  J  (£u  +  mv)p  ds 


(5) 


whi  ch  applies  to  the  body  of  liquid  enclosed  in  a  cylinder  of  unit  length 
In  the  z  direction  and  having  s  as  the  circumferential  boundary  curve. 
Here  i  and  m  are  the  directional  cosines  of  an  inwardly  directed  normal 
to  the  bounding  curve  element  ds  or  the  bounding  surface  of  area 
dA  =  1  *  s  .  The  quantities  T  and  V  refer  to  kinetic  and  potential 
energies  of  the  liquid  bounded  'by  s  and  are  determined  from 


and 


T  = 

V  = 


u2  +  v2)dS  ;  dS  =  dx 
dS  ;  dS  =  dx  dy 


dy 


Multiplying  the  two  sides  of  equation  5  by  dt  and  integrating  between  two 
instants  t.  and  t2  yields 

(T  +  V)„  -  (T  +  V)..  =  I  f  (iu  +  mv)p  ds  dt  (6) 

J\Js 

The  interpretation  is  that  the  total  increase  in  energy,  potential  and 
kinetic,  of  the  same  portion  of  liquid  in  its  path  of  displacement  is  equal 
to  the  work  done  by  the  pressures  on  its  surface.  This  is  stated  by  Lamb. 
See  fig.  2a.  Since  the  expression  in  equation  6  gives  the  change  in  the 
energies  of  the  same  portion  of  liquid  in  its  path  of  displacement  for  the 
times  t1  and  t2  ,  this  expression  may  be  referred  to  as  the  Lagrangian 
x’orm  of  the  equation  of  energy. 

12.  An  alternate  form  of  the  equation  of  energy  is 


where 


lo  •  L ty  and  3 
unaf  the  change 


v  and  s  represents  a  curve  drawn  in  the  field  of  ve- 
is  the  area  enclosed  by  the  curve.  The  interpretation  is 
in  kinetic  energy  in  a  region  delineated  in  tne  flow  field 
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by  a  closed  curve  is  due  to  the  inflow  and  outflow  or  potential  energy  and 

k i  n^t.i  f  faripv/ru*  intrt  wnn  nunv  fmm  t.h^  enclosed  ?.rO?.j  and  t-0  the  WOrK 

dona  by  the  pressures  on  the  body  of  liquid  in  the  delineated  area.  See 
fig.  2b.  Since  the  expression  in  equation  7  gives  the  changes  in  the 
energies  in  the  fixed  area,  this  expression  may  be  referred  to  as  the 
Eulerisn  form  of  the  equation  of  energy.  The  proof  of  the  expression  in 
equation  7  will  be  given  in  the  Addendum. 

13.  In  the  ordinary  application  of  the  energy  relations,  usually  the 
flow  is  through  passages  with  rigid  walls,  for  which  at  the  walls  u  and 

v  vanish. 

Pneumatic  Pressures 

14.  Use  may  now  be  made  of  equation  6  to  examine  the  increase  in 
wave  energy  with  time  in  the  channel  downstream  of  the  nozzle  mouth.  See 
fig.  lb.  It  is  assumed  that  at  time  t  <  0  the  liquid  in  the  channel  is 
at  rest;  the  wave  motion  commences  at  time  t  =  0  ;  and  during  time  t  =  t 
the  front,  a  node  point,  has  traversed  the  distance  x  =  xQ  . 

As 

kXQ  =  <rt  (3) 

the  disposition  of  the  surface  of  the  wave  is  . 

h  =  h  sin  (kx„  -  kx)  ;  0  <  x  <  xn 
mu  u 

and 

h  =  0  ;  x  >  xQ  (9) 

15.  Consider  the  area  under  the  wave  bounded  by  the  wave  surface, 

the  vertical  A’B’  at  x  =  0  ,  the  vertical  CD  at  x  =  Xq  ,  and  the  bot¬ 
tom  segment  B'D  .  Call  this  the  final  area  .  The  liquid  of  this  area 
initially,  that  is  at  time  t  =  0  ,  was  contained  in  the.  area  of  the  rec¬ 
tangle  bounded  by  the  vertical  AB  at  x  =  -Ax  ,  the  straight  line  AC  , 
the  vertical  CD  at  x  =  x^  ,  and  the  bottom  segment  BD  .  Call  this  the 
initial  area  S^  .  Since  the  liquid  is  incompressible,  and  are 

equal  and 
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h  dx 


8 


(10) 


and  S.S  kx^  ~  at  s  also 


pgh  IL  Dgh 

Eg  =  0'^'  (cos  at  -  l)  *’•■  ■■g-”  (at  -  sin  at  cos  at) 


l6.  This  energy  comes  from  the  work  dene  by  the  pressure  ,  pre¬ 
vailing  at  the  section  x  =  0  ,  on  the  line  AB  ,  and  the  work  is 

Ei  -fofo°  P°‘  dy  dt 

The  integral  will  equal  Eg  only  if 

PQ  =  pgh  +  pgHQ  -  pgy  =  pg(HQ  +  h  -  y)  (il) 


where  h  is  the  surface  displacement  at  x  =  0  .  Recalling  that 

u  =  u  sin  at  and  h  =  h  sin  at  .  the  integral  for  E,  yields 
m  m  i 


E, 


u  dt 

=  (cos  at  -  1)  +  (at  -  sin  at  cos  at-) 


As  u^kHg  -  ahm  ,  see  equation  1,  also 

2 

Pgh  H.  pgh 

Ei  =  — g~  (cos  at  -  1)  +  (at  -  sin  at  cos  at)  (12) 

The  equality  between  E^  and  Eg  ,  thus,  is  shown. 

17.  One  recognizes,  of  course,  that  the  pressure  PQ  prevailing  at 
the  nozzle  mouth  (that  is,  at  section  x  =  0)  is  made  up  of  two  parts,  one 
due  to  the  wave  disturbance  h  and  the  other  due  to  the  undist’irbed  liquid 
o£  depth  Hq  .  Having  ascertained  the  magnitude  of  PQ  ,  equation  11,  the 
proper  pressure  P  in  the  pneumatic  chamber  over  the  water  surface  may 
next  be  evaluated  after  resorting  to  the  alternate,  or  the  Eulerian,  form 
of  the  equation  of  energy,  equation  7«  For  the  present  case  (see  fig.  lc), 
the  bounding  curve  S  io  made  of  the  following  parts:  the  segment 
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*onsisiIng  of  the  free  surf;  o  of  the  liquid  in  the  chamber,  the  segment 
S0  to  be  identified  as  the  nozzle  mouth,  and  the  rigid  boundaries  of  the 
chamber  and  the  nozzle  in  contact  with  water.  The  appropriate  form  of  the 
energy  equation,  then,  is 


+  PQ)dy  (7a) 


where  denotes  the  water  surface  velocity  in  the  chamber  and  Up  the 

velocity  at  the  mouth.  The  sums  of  the  terms  in  the  parentheses  of  the 
individual  integral.,  are  constant,  and  since 


V  =  v*o 

the  simplified  fonr-  of  the  energy  equation  is 

JJ 

fa  5t“  d:<  dy*S  Vdy  =  (2  vi  +  ^1  +  p)  ”  (2  *2  +  ^2  +  po) 

V 

18.  For  the  purpose  of  evaluating  the  term  on  the  left-hand  side, 
measure  first  the  vertical  distances  in  term^  of  ,  and.  the  horizontal 
distances  in  the  direction  of  wave  motion  in  terms  of  \/2k  .  Thus, 
introducing 


x'  =  x/k  and  y'  -  2«y/X  =  y/k 
the  velocity  q  now  is  expressed  as  the  function 

q  =  u  f(x' ,y‘)  sin  at 

and  with  it  the  left-hand  term  in  equation  7b  may  be  evaluated,  giving 
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■IW' J. ■■urium »  ujli  *” 8T*T 


«£ 


jTf^  “  '*  :fg  V*  -  »t  2  f  u»  cos  ot  “  Ht  f  s»m  “s  ot 


U3) 


where  II,  is  a  pure  nunber,  its  value  depending  on  the  shape  of  the  nozzle. 

Xf 

Using  this  result  and  next  expressing  in  terns  of  Ug  =  (u^  sin  •&), 

the  energy  relation  equation  7b  simplifies  to 

2 


2  jNtghm  cos  ot  +  u2  1  -  (-^)  sin2  ot]J  =  nf^  -  oft. 


+  P  -  P 
2  r  0 


Introducing  into  the  right-hand  member  the  va?ue  of  Pq  from  equation  11, 


\  -  s(^q  +  S)  ,  n2  =  gy 


and  also  writing 


P  =  pgh^  -  PgS 


m 


the  final  form  of  the  energy  becomes 


m  +  \  VHa  cos  +  I  \  .!  *  (t)  ]  Si 


sin2  at 


Since  h  is  small  in  comparison  with  H„  ,  the  last  term  on  the  right- 
flu  2  ^ 

hand  side  involving  will  be  small  in  comparison  with  gh  and,  thus, 
will  be  ignored.  One  now  has 


and  according!;  •  the  desired  pneumatic  chamber  pressure  should  be 


—  =  h  sin  ot 
pg  m 


(  Kt  ^ 

-  —  V  cos 


(15) 


This  may  be  put  n  the  form 


—  =  A  sin  (ot.  -  oe) 
Pg  '  ' 


(16) 
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£ 


f 


K- 


& 


E 


I 


where 


and 


the  value  c:'  H.  being  evaluated  from  equation  13. 

v 

19.  One  may  now  summarize  the  results.  The  generation  of  long  waves 
of  wave  height  2;^  and  of  period  T  will  be  brought  about  in  a  practical 
manner  through  the  introduction  of  sinusoidally  varying  pressure  P  in  the 
pneumatic  chamber  of  period  T  and  of  semi  amplitude  pgf  4  The  phase  dif¬ 
ference  in  the  pressure  variations  inside  and  the-  wave  surface  oscillations 
outside  is  given  by  e  .  Pressure  P  denotes  the  excess  or  deficiency  in 
reference  to  atmospheric  pressure. 

20.  In  the  analysis  given  above  and  leading  to  the  chamber  pressures, 
losses  in  th  chamber  and  nozzle  passages  are  ignored.  The  theoretical 
treatment  of  the  losses  would  certainly  be  sanevhat  uncertain  for  nozzles  of 
arbitrary  shape.  To  account  for  the  losses,  one  may  increase  the  semi  ampli¬ 
tude  of  the  pressure  variations  A  by  an  amount  5A  .  If,  in  the  nozzle 
used,  the  two  sides  of  the  throat  are  not  s^n metrical,  losses  would  be  of 
different  amounts  for  the  inflow  into  the  chamber  and  the  outflow.  As  a  re¬ 
sult.  waves  generated  in  the  proximity  of  the  chamber  would  contain  higher 
harmonics.  Perhaps  this  difficulty  can  be  overcane  by  choosing  a  nozzle 
having  small  3. esses.  Ibis  is  a  matter  that  could  be  examined  in  models. 

21.  Simultaneously  with  the  development  of  thought  and  the  design  of 

pneumatic  wave  genera  -s  at  California  Institute  f  Technology,  similar 

developments  were  c  sd  at  the  David  Taylor  Model  Basin  leading  to 

4 

elaborate  machines  for  the  generation  of  deep-water  waves.  In  the  latter 
type  machines  the  diffuser  nozzles  are  avoided.  Under  the  direction  of 
Mr.  Brownell,  a  system  of  generators  with  accessory  controls  has  been 
evolved  which  is  capable  of  simulating  a  composite  sea  state  of  any  speci¬ 
fied  complexity.  The  design  criteria  of  the  machines  were  established  frem 
models . 
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22.  In  these  types  of  machines  with  nozzles  absent,  if  desired,  the 
mechanics  of  flow  may  be  adequately  described  with  results  from  the  momen¬ 
tum  principle.  A  form  of  the  analysis  along  these  lines  is  given  by 

5 

Kergoat  in  a  brief  communication.  In  the  paper  are  included  experimental 
results  which  show  good  agreement  with  theory. 
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PART  m:  APPROXIMATE  THEORY-  OP  IKK3IA1IC  GEKERAT0R3  FOR  SURGES 

23.  The  feasibility  of  pneumatic  wave  generators  for  the  purpose  of 
creating  a  singular  positive  translation  wave  of  great  length  was  being 
examined  experimentally.  In  this  section  an  approximate  theory  of  the 
generator  is  developed  to  serve  as  a  guide  for  the  analysis  of  the  test 
data. 

Wave  Characteristics 

24.  It  is  desired  to  produce  an  elongated  wave  of  constant  height  h 
advancing  with  a  constant  velocity  of  propagation  ix.  still  water  of 
constant  depth  Hq  .  Applying  the  Boussinesq  concept  of  propagation  to  the 
present  case 


ox 


(17) 


sirce  d2h/dx2  vanishes. ^  Let  Ug  be  the  particle  velocity  in  the  wave, 
constant  in  any  vertical  section.  By  the  continuity  condition 


UgO^j  +  h)  »  oh 

and  hence,  neglecting  powers  of  h/tiL  , 

Eliminating  to  ,  using  equation  17, 

Operation  of  Wave  Generator 


(18) 


(19) 


2$.  The  generator,  which  consists  of  an  air  chamber  with  a  long  dif¬ 
fuser  type  nozzle  at  the  base  opening  into  a  channel,  is  shown  in  fig.  3*» 
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This  is  of  the  type  to  be  referred  to  as  a  low  generator.  At  the  top  are 
two  apertures,  one  of  them  connected  to  an  aspirator  and  the  other  to  out¬ 
side  air.  The  openings  are  controlled.  The  pressure  in  the  pneumatic  cham¬ 
ber  is  changed  from  pQ  to  pQ  -  £spQ  ,  where  £f>Q  is  a  positive  quantity, 
by  opening  the  aperture  leading  to  the  aspirator.  As  water  in  the  chamber 
is  raised  to  a  level  H  ,  and  water  remaining  in  channel  is  of  depth  Hq  , 

P«H  -  ^5>0  -  cgHQ 
or 

%  =  Pg(K  -  H0)  (20) 

At  t  =  0  the  aspirator  is  disconnected,  and  the  aperture  to  the  outside 
is  opened,  allowing  air  to  enter.  At  time  t  let  the  chamber  relative 
pressure  be  *  6  the  fall  of  water  surface  in  the  chamber,  L  the 
length  of  the  wave  generated.,  and  h  the  height  of  the  wave  with  respect 
the  undisturbed  water  level. 

26.  If  is  the  downward  velocity  of  the  falling  surface, 

v1  *  d6/dt  (21) 

By  the  condition  of  continuity 


where  1  is  the  longitudinal  length  of  the  air  chaaber.  Differentiating 
with  respect  to  t  and  since  dL/dt  =  © 

VjJ  *  <uh 

Also 

Tl*  *  *2^) "  h) 

Problem  of  Chamber  Pres  sure 

27.  Proa  the  above  developments  it  is  seen  that  a  requisite  for  the 


(23) 

(2b) 
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generation  of  a  wave  of  constant  height  h  moving  with  a  constant  speed 
of  propagation  oj  in  a  channel  of  constant  depth  Hq  is  that  the  water 
surface  in  the  pneumatic  chamber  must  fall  uniformly  (that  is,  with  a  con¬ 
stant  rate).  This  poses  a  restriction  on  the  reduced  pressure  £p  ,  and 
the  dependence  of  on  t  may  be  determined  conveniently  by  again 

using  the  principle  of  energy. 

28.  In  fig.  3a  the  disposition  of  liquid  in  the  channel  to  the  left 
of  the  limiting  section  A-A'  ,  and  in  the  pneumatic  chamber  for  the  ini¬ 
tial  time  t  =  0  ,  is  shown.  In  fig.  3b  the  subsequent  disposition  of  the 
same  liquid  at  time  t  is  shown.  The  limiting  section  A-A'  is  at  such  a 
distance  from  the  mouth  that  the  issuing  wave  has  not  reached  it  at  the 
•*  ime  t  .  Consider  the  energies  of  the  water  to  the  left  of  A-A*  for 
these  two  instants,  t  =  0  and  t  =  t  .  Since  at  time  t  =  0  the  liquid 
is  at  rest,  the  initial  energy  cf  the  system  is  only  potential.  On  the 
basis  of  the  notations  of  fig.  3a, 

(T  +  V)x  =  (H2*  +  (25) 

At  time  t  the  liquid  has  kinetic  energy  made  of  two  parts,  one  (T21)  re~ 
lating  to  the  water  in  the  chamber  and  in  the  nozzle,  and  the  other  (Tpp) 
.•elating  to  the  wave  in  the  channel.  Computation  will  show  that 

*21  *  I  A  (» - 8  -  V*  ♦  1 (26) 


where  is  a  numerical  constant  the  value  of  which  depends  on  the  shape 

of  the  nozzle  end  the  path  connecting  the  min  body  of  water  in  the  chamber 


with  the  nozzle.  Between  the  levels  of  H  and  Hg  the  cross  section  of 
the  pneumatic  chamber  remains  the  same.  With  low  generators  as  in  fig.  3, 


Ilg  is  identical  with  the  narrow  opening  D  of  the  nozzle.  The  expression 
of  N+  in  terms  of  nozzle  dimensions  is  derived  in  the  Addendum  for  two 


types  of  pneumatic  generators,  the  low  and  the  elevated.  Ac  this  time  the 
exact  numerical  value  of  is  not  needed.  As  regards  the  second  part 

of  the  kinetic  energy,  one  has 
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T22  =  |  ^  (H0  +  h)L 


(27) 


where  L  is  the  length  of  the  wave.  Finally,  the  total  kinetic  energy, 

T2  =  T21  +  T22  *  is 

Ta  ‘  f  \  (»  *  6  -  V  +  f  “^t  +  I4H0  (X  +  ^)  (28> 

29.  The  potential  energy  of  the  liquid  to  the  left  of  A-A’  in  the 
chamber  and  in  the  channel  for  the  instant  t  is 


V2  =  f  (K  -  &)2i  +  f-  (Hq  +  h)2L  +  &  (Lw  -  L)H^ 


(29) 


Write  (T  +  V)2  for  T,-,  +  Vp  .  Denote  the  difference  of  the  energies  for 
the  instants  t  =  t  and  t  =  0  as 


ZJE  =  (T  +  V)2  -  (T  +  V)., 


(30) 


Substituting  from  the  above  and  also  making  use  of  the  continuity  condition 
6/  =  Lh  ,  since  h  is  constant  throughout, 

&  =  -Cg(K  -  Hq)M*  +  f  «6hl  +  §  (gh2  +  U^L 

+U4- \)hL + f 4 <H -  y* +  i 4&t  <31) 

2 

For  a  further  simplification  of  this  expression  one  may  first  neglect  v, 

2 

in  comparison  with  Ug  ,  since  1  is  many  times  larger  than  .  Again, 
two  sets  of  terms  may  be  grouped  separately  as  follows 


|(u|h  +  gh2  +  u^)  =  gh2N 


(32) 


where 


h“1  +  H 
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and 

\  (K  -  Hgji  +  u^Nt  =  ghl§4  (33) 

where 


30.  Tne  simplified  expression  for  £E  ,  equation  31>  now  is 

AE  =  pg[|  hH^M  +  tfa  -  (H  -  KqJKL  +  |  6MjJ  (3*0 

Since  L  =  ut  and  5  =  v^t  ,  the  last  may  be  written  also  as 

Z5E  -  pg[!  hH^l  +  h2^  -  (H  -  HQ)hL  +  ~  6M,]  (35) 

Here  ZSI  represents  the  difference  between  the  energy  of  the  system  at 
time  t  and  that  at  time  t  =  0  .  The  quantities  M  and  N  are  inde¬ 
pendent  of  time  and  have  set  values  for  H  ,  Hq  ,  and  h  constant;  and, 
thus,  dE  is  a  quadratic  algebraic  function  of  time.  The  variation  of  Ap 
with  time  must  be  such  as  to  account  for  this  particular  form  of  . 

What  the  desired  variation  is,  as  suggested  before,  may  be  derived  from  the 
energy  relation  in  equation  6.  If  W  is  the  work  done  by  the  pressures  on 
the  boundaries  of  the  liquid  to  the  left  of  A-A’  ,  by  the  energy  principle 


where 


££  =  W 


tM  dt 


(36) 

(37) 


Since  v^ 


is  constant  along  the  water  surface  in  the  pneumatic  chamber 


W  = 


4 

~  r\ 


v.j£p  dt 


(38) 


Before  proceeding  further  let  us  consider  the  physical  conditions  at  the 
instant  of  wave  generation.  At  t  =  0  ,  the  water  in  the  pneumatic  chamber 
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is  at  rest.  At  t  =  e  and  thereafter,  water  particles  in  the  chamber 
will  be  moving  with  constant  velocity  v^  ,  since  it  is  assumed  that  for 
t  >  e  the  wave  is  of  constant  height.  On  this  basis 

*€ 


W  =  -lfo  Y^*  it-  tv^  Cf 


dt 


(39) 


Here  v.  1  denotes  the  surface  velocity  of  water  in  the  chamber  for  times 

Xj. 

less  than  €  .  We  may  write  also 

re  c1 

W  a  -iv^  I  -XX- - i  Ap  dt  -  iv ^  I  £ p  dt 

•'o  1 

or 


and  if  we  put 


W  -  -iv^0(-lc  *T  jg-it) 


(bo) 


(41) 


Using  one  of  the  forms  of  the  condition  of  continuity,  equation  2k,  and 
introducing  the  value  of  £pQ  ,  equation  20,  yields 

W  =  -u^  ♦  h)p.S(H  -  Hg)  (-le  +J  ^  dt  j 


or 


W  =  -p0ch(H  -  1^)  ^-l€  + J  dt^ 

31.  For  t  >  e  ,  we  suppose  that  the  pressures  sure  given  as 


(42) 
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T&-  ~  1  -  *i  *  *2* 


and  if  we  denote  by  the  pressure  that  would  exist  at  t  =  0 


(43) 


and  hence 


How 


and  this  makes 


^i  =  U  - 


dt  =  (l  -  k1)t  -  kgt2 


m 

(45) 


W  =  -pgmh(K  -  Hq)  [-Ie  +  (l  -  k^t  -  l^t2]  (46) 

which  like  £E  ,  equation  35-,  is  an  algebraic  quadratic  function  of  t  . 
Comparing  terms  in  the  two  expressions,  equations  35  and  46,  not  involving 
t  one  finds 


pg(H  -  H^Ie  =  |  (47) 

Accordingly,  since  M  is  a  positive  number,  I  is  also  a  positive  number. 
The  suggestion  is  that  there  should  be  a  relatively  sudden  increase  of 
pressure  which  then  after  a  short  period  of  time  falls  to  .  Further 

discussion  as  regards  the  transient  pressure  may  be  omitted  here.  Initial 
conditions  will  be  discussed  later. 

32.  Comparing  the  coefficients  of  t  in  the  expressions  of  and 
W  ,  equations  35  and  46,  one  finds 


or 


hN  -  (H  -  Hq)  =  -(1  -  kjXH  -  Hq) 


(43) 


«  (H  -  Hq)!^ 
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This  expression  enables  one  to  determine  in  terms  of  the  wave  height 

h  and  the  elevation  of  the  water  surface  in  the  pneumatic  chamber  H  at 

the  time  the  wave  starts.  From  one  next  obtains  the  value  of  the 

pressure  in  the  chamber  to  generate  the  wave  of  the  desired  height  h  . 

2 

33*  Comparing  the  coefficients  of  t  in  the  expressions  of  Z® 
and  W  ,  equations  35  and  46,  one  finds  that 


=  (H  -  H^kg 


(50) 


From  the  continuity  relation 


Th~  above  yields 


which  determines  what  the  rate  of  fall  of  pressure  in  the  pneumatic  chamber 
should  be  in  order  to  maintain  a  constant  wave  height  with  time  in  tile  gen¬ 
erated  wave.  Here  kg  is  determined  in  terms  of  h  ,  H  ,  and  £  ,  the 
latter  being  the  dimension  of  the  tank  in  the  longitudinal  direction. 

34.  Equation  50  may  be  written  in  another  form.  Since  v^  =  d&/dt 
and  2kg  =  -d(ZJp/fip0)/dt  , 

i  =  <H  -  V  £(%) 

But 


£p0  =  pg(H  -  Hq) 


and  hence 
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(50a) 


dfr  d  /£%>  \ 
dt  ~  "  dt  V  pg  ) 

indicating  that  the  rate  of  fall  of  water  surface  in  the  chamber  equals  the 
rate  of  increase  of  chamber  pressure. 

35.  The  relation,  equation  49  5  connecting  h  with  is  for  the 
condition  that  in  the  flow  of  '.niter  through  the  nozzle  there  is  no  loss  of 
energy  from  friction  associated  with  turbulence  and  deflected  motions. 

Since  in  the  best  designed  nozzle  some  loss  is  expected  to  occur,  a  provi¬ 
sion  must  be  made  to  account  for  the  effect  of  the  loss.  The  frictional 
less  reduces  the  energy  available  to  create  the  wave.  Let  this  loss  be 
£L  ,  and  write 


=  Xpgh^HL 


(52) 


as  discussed  in  the  Addendum.  This  means  that  in  the  energy  expression 
fsE  ,  equation  37,  in  the  place  of  pgh^RL  one  must  now  write  (l  +  X)pgh?ML . 
Proceeding  as  before  one  now  has 


h_  __  /  _H 


(53) 


We  may  refer  to  X  as  the  nozzle  friction  facto  . 

36.  Examining  the  two  expressions,  equations  51  and  53»  it  is  seen 
that  the  quantities  k^  and  kg  i— e  not  independent  but  are  related. 
Comparing  these  equations  with  each  other,  it  is  found  that 


2k2  _  _1__  l  h\ 

^"1  +  X  i  \  +  2  Hq/ 


(54) 


so  that,  if  the  friction  factor  is  known,  ohe  ratio  2kg/k^  can  be  evalu¬ 
ated  from  the  assumed  values  of  h  and  Hq  and  the  chamber  size  1  . 
Conversely,  if  2kg  and  are  observed,  the  relation  is  utilized  to 
compute  X  ,  since 


1  +  X 


h'T&o/. 

2*2  ‘  V 


(55) 
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Experimental  Evidence 

37.  The  preceding  analysis  dealing  with  the  generation  of  an  elon¬ 
gated  wave  of  constant  height  suggests  that  this  kind  of  wave  would  he  more 
readily  and  easily  generated  if  the  pneumatic  chamber  is  elevated.  If  the 
chamber  is  elevated,  the  total  fall  of  water  surface  therein  would  be  a 
smaller  fraction  of  the  initial  height  H  ,  and  accordingly  the  pressures 
could  be  conveniently  prodc/.  ed  without  a  variable  control  of  the  air  open¬ 
ing.  To  examine  the  feasibility  of  the  idea,  experiments  were  conducted 
with  a  generator  as  shown  in  fig.  4.  The  relations  shown  in  equations  51 
and  53  still  hold.  The  only  modification  would  be  in  regard  to  the  coef¬ 
ficient  M  ,  equation  47,  which  changes  with  the  type  of  generator  adopted. 
This  is  not  important  since  transient  pressure  effects  at  the  time  of  wave 
issuance  will  be  ignored. 

33.  For  the  present,  the  important  observation  data  to  be  examined 
are:  the  wave  height  at  the  mouth  of  the  nozzle  as  a  function  of  time,  the 
pressure  variation  in  tie  chamber,  and  the  fall  of  water  surface  inside  the 
chamber.  An  example  of  a  tracing  from  an  electronic  record  is  shown  in 
fig.  5.  The  lower  curv?  shows  the  variation  in  wave  height  in  the  vicinity 
of  the  noznle  mouth.  Note  that  for  a  long  time  the  wave  height  remains 
constant.  The  middle  curve  shows  the  increase  with  time  of  the  pressure  in 
the  chamber  and  in  the  area  above  the  water.  If  one  ignores  the  oscilla¬ 
tion  of  pressure  in  the  short  dura l ion  following  the  instant  that  the  air 
passage  is  opened,  the  increase  of  the  pressure  with  time  is  linear.  It 
would  be  adequate  to  represent  the  pressure  course  by  a  single  straight 
line  for  all  instants.  Note  that  the  intersection  of  tne  pressure  line 
with  the  zero  time  axis  is  indicated  by  .  For  the  present  the  reason 

for  the  pressure  oscillations  for  small  t  values  will  not  be  dis  ussed. 
The  matter  will  be  taken  up  in  another  section.  The  upper  curve  shows  the 
fall  of  the  water  surface  in  the  chamber  with  time.  The  rate  of  fall  is 
constant,  and  its  absolute  value  equals  the  rate  of  rise  of  the  pneumatic 
pressure.  These  results  are  in  accordance  with  theory.  Tne  data  obtained 
from  this  and  a  few  other  records  are  presented  in  table  1. 

39*  The  data  will  be  examined  to  compar  »  the  theoretical  results 
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with  the  observed  results.  The  quantities  involved  are 


which  is  the  expression  for  the  particle  velocity  in  the  wave  ; 


(19  bis) 


(51  bis) 


which  is  the  expression  for  the  rate  of  increase  of  pressure  in  the  pneu¬ 
matic  tank  and  time;  and 


(53  bis) 


which  is  the  expression  that  determines  the  required  (effective!)  pressure 
in  the  tank  at  the  time  when  the  wave  commences  to  issue  from  the  nozzle. 

ho.  The  constants  and  X  appear  in  the  last  equation,  and  thus 
it  would  serve  to  determine  X  ,  the  nozzle  resistance  factor,  in  terms  of 


h  and  as  observed  quantities.  Another  way  to  determine  X  is  from 


.  .  .  _  *1  ^0  /.  .1  h\ 
1  +  x  =  au  t  V  2  L  / 


(55  bis) 


4l.  It  is  helpful  to  mention  that  the  quantities  k^  and  kg  are 
not  observed  but  are  derived  from  the  obse  ved  quantities  Ap^/pg  ,  Ap/pg  , 
and  d/dt  ♦  Ap/pg  ,  using  the  definitions 


and 


2*2  -  -  ^  tofeo 

k^  =  1  -  Apj/Apg  (See  equations  h.3  and  45.) 


42.  The  computed  values  are  shown  in  the  lower  part  of  table  1. 
There  is  a  fair  degree  of  agreement  between  the  observed  and  the  computed 
values  of  Ug  and  kg  .  There  are  two  determinations  of  X  ,  and  as  ex¬ 
pected,  there  are  considerable  variations  from  one  run  to  another. 
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PART  IV:  THE  MECHANICS  OF  PRESS' TRE  BUILDUP  IN  PNEUMATIC  CHAMBERS 


43.  It  was  shown  in  Part  IH  that  the  mere  raising  of  the  level  of 
the  water  in  the  chamber  and  allowing  the  outside  air  to  enter  through  an 
aperture  without  the  assistance  of  additional  controls  is  sufficient  to 
cause  the  pressure  to  fall  at  a  constant  rate.  The  fall  is  as  required  by 
the  theory  of  long  waves  of  constant  height.  This  behavior  of  the  gener¬ 
ator  needs  to  be  explained  on  the  basis  of  the  mechanics  of  air  flow  into 
the  chamber. 


Mass  Rate  of  Air  Inflow 


44.  As  a  first  step  in  the  analysis,  consider  the  closed  vessel 
shown  in  fig.  6  which  has  an  opening  in  the  form  of  a  circular  orifice  of 
area  a  .  Let  the  volume  of  the  vessel  be  Vq  .  When  the  pressure  inside 
is  below  atmospheric  there  is  a  mass  rate  of  flew  inward.  Let  the  condi¬ 
tions  outside  be  given  by  pQ  and  pQ  and  those  of  the  interior  (that,  is, 
at  the  vena  contracts)  p^  and  .  Let  be  the  area  of  the  vena  con¬ 
tracts.  Assume  an  adiabatic  flow.  For  the  mass  rate  of  flow  one  has  from 


where  is  the  velocity  of  air  at  the  vena  contracts,  is  the  veloc¬ 
ity  of  sound  in  the  air  outside,  that  is 

co  *  (7  *  (»«o>1/6  t57> 

and  7  is  the  ratio  of  the  specific  heats  at  constant  pressure  and  con¬ 
stant  volume.  For  air  7  =  1.408  ,  and  for  standard  atmosphere 
Cq  -  IO89  ft/se*  ;  is  the  temperature  in  degrees  Rankine  corresponding 
to  32  F. 

45.  In  our  applications  the  difference  between  pQ  and  p^  is 
small.  Denoting  the  difference  by  dp  »  Pj/Pq  *  1  *  ^?/Pq  • 
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One  nay  refer  to  C  as  the  coefficient  of  discharge,  since  it  is  a  dimen¬ 
sionless  quantity. 

46.  To  coefficient  C  one  nay  ascribe  a  meaning  more  general  than 
that  implied  by  equation  58,  Equation  59  "By  be  used  also  in  the  cases 
where  there  is  dissipation  in  the  flow  reaching  the  vena  contracts  or  the 
passage  may  be  through  a  tube  of  any  kind.  The  formula  serves  to  relate 
the  mass  flow  of  air  to  the  pressure  in  the  closed  chamber.  In  such  cases 
C  may  be  thought  to  depend  on  Reynolds  number  Be  =  V?  •  ^  ,  d  being  the 
diameter  of  the  tube. 

47.  It  is  advantageous  to  express  the  relation  in  equation  59  in  a 
dimensionless  form.  Omitting  the  steps  of  transformation  the  expression  is 


(5>9«0 


To  illustrate  we  now  consider  the  Fliegner  formula  cited  by  Praadtl. 
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Poison  et  al.  have  shown  the  validity  of  the  formula.  In  metric  units 

Q  =  0.76a  -J  (px  - 

In  the  terms  of  this  paper  and  if  £sp  is  small  with  respect  to  pQ  , 

g  H  =  0.76a  .  Pq/6 

One  has,  omitting  the  steps  of  transformation, 

dm  „  ^  R1/2  .  /Ap  .  C0aC0 .  /Ap 

=  0.76a  — pa  \lxr  =  1.31  — T7o\  ?r 
dt  6  0  V  P0  v1/2  V  P0 

since  R/g  =  29.4  meters  per  degree.  Accordingly,  C  =  1.31  and  the  vena 
contracts  occupies  the  entire  area  of  the  opening. 

Determination  of  Coefficient  of  Discharge 

48.  In  the  applications  to  be  made  eventually,  it  will  be  necessary 
to  know  the  value  of  the  discharge  coefficient.  A  method  to  determine  it 
would  be  to  ascertain  the  rate  of  increase  of  the  pressure  inside  the  ves¬ 
sel  with  an  initial  difference  of  the  pressure  ApQ  .  If  the  volume  VQ 

is  large  enough  the  filling  may  be  regarded  as  an  adiabatic  process.  Then, 
the  pressure  and  density  in  the  vessel  are  those  at  the  vena  contracts. 

49.  In  this  hypothesis  one  uses 


and 


m  =  0lV0 


35 i  =  v 

dt  dt  0 


(60) 


since  Vq  remains  constant.  Now  for  adiabatic  conditions 


b  ifp  /piV 

%  V°o/ 


(61) 
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Sin-re  Ap  is  a  small  quantity 


1 

7p0 


and  differentiating  with  respect  to  time 


fjO  ffl  _  _  1  dZ£ 
oQ  dt  "7  dt 

Introducing  the  sound  velocity  cq  >  cq  =  TPq/Pq 

ffl  _  1  dAp 

dt  "  ”  2  dt 
C0 

and  from  equation  60 

V 

dm  0  dAp 
dt  "  2  dt 

“0 

Equating  this  to  the  mass  rate  of  entry  from  outside,  equation  59, 


or 


cr 


0  d  A 
-dt^  = 


Q  Ap 

dt  Ap0 


(62) 


2S 


The  solution  is 


50.  An  experimental  determination  of  C  was  carried  out  by  noting 
the  filling  rates  of  a  rectangular  box  2  ft  high  and  with  a  cross  s^-Uon 
of  1  by  2  ft.  In  the  upper  plate  (0.485  in.  thick)  numerous  circular  aper¬ 
tures  with  varying  diameters  were  drilled.  The  edgc-s  of  tie-  apertures  w----r-- 
sharp.  Taking  an  aperture  at  a  time  when  the  air  pressure  in  the  vessel 
was  reduced  by  an  amount  Ap^  ,  the  aperture  was  next  opened  to  the  air 
outside,  and  the  increase  of  the  pressure  within  the  vessel  was  noted. 

Fig.  7  is  an  example  of  a  plot  showing  the  linear  variation  of  (  Ap/Ap^ )  ! 
with  time.  This  is  in  agreement  with  the  theoretical  form  in  equation  63. 
Thus,  if  m  is  the  slope,  the  coefficient  is 

C  =  2n0Ca 

The  values  of  C  for  orifices  of  various  diameters  are  shown  in  table  2. 
The  coefficient  of  discharge  is  practically  independent  of  the  diameters 
except  that  the  coefficient  shows  a  slight  increase  in  value  when  the  ratio 
of  plate  thickness  to  orifice  diameter  is  increased.  The  constant  value  is 
C  =  0.71. 

51.  The  tests  were  continued  with  vessel  volumes  VQ  smaller  than 

4.03  ft.  These  lesser  volumes  were  conveniently  obtained  by  filling  the 

original  vessel  with  water.  The  results  of  the  tests  are  shown  in  table  3. 

It  is  seen  that  C  when  determined  by  equation  63  decreases  with  Vr  .  In 

o 

the  routine  tests  on  wave  generation  it  was  observed  that  owing  to  inertia 
the-  water  surface  in  the  chamber  remains  still  for  a  short  time  following 
the  instant  when  the  chamber  is  opened  to  outside  air.  This  would  mean 
that  for  the  initial  period  the  air  is  flowing  into  a  closed  vessel  of  un¬ 
varying  volume.  In  figs.  8  and  9  the  initial  increase  of  pressure  with 
time  for  the  two  types  of  generators  is  shown.  The-  air  openings  were  in 
the  form  of  circular  orifices.  Since  the  manner  of  the  increase  is  in 
agreement  with  equation  63,  C  may  be  computed  by  this  formula.  The  * 
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ro cults  of  the  reductions  are  shown  in  table  4. 

52.  It  will  be  shown  in  the  Addendum  that  the  effective  values  of 
C  (that  is,  the  values  computed  on  the  basis  of  equation  63)  may  be  depen¬ 
dent  on  the  volume  V.  occupied  by  the  air  jet  and  the  length  X.  of  the 

0  J 

jet.  provided  that  a  short  distance  below  the  air  opening  the  jet  is  de¬ 
flected  laterally  over  the  surface  of  the  water.  This  would  suggest  that 
C  could  be  a  function  of  V./V0  and  X./A  where  A  is  the  initial  air 
gap  between  the  water  surface  and  the  chamber  ceiling.  In  fig.  10  the 
effective  discharge  coefficient  is  plotted  against  X./A  .  There  is  a  fair 
congruency  of  the  points  with  the  curve  drawn.  For  the  cases  considered, 
the  effect  of  the  ratio  V./V0 


on  C  is  negligible.  The  meanings  of  X. 

J 


and  V.  are  discussed  in  the  Addendum. 

J 

53.  In  a  more  general  sense,  as  previously  mentioned,  C  will  be 
interpreted  to  be  a  proportionate  factor  in  the  empirical  formula,  equation 
59)  connecting  the  mass  rate  of  flow  of  air  with  the  pressure  in  a  receiv¬ 
ing  chamber.  The  air  passage  may  be  of  any  form.  In  the  first  runs  of 
pneumatic  wave  generation  the  chambers  were  connected  to  the  aspirator  and 
to  the  ail  outside  by  means  of  a  so-called  three-way  valve.  In  fig.  11  are 
shown  the  initial  rises  of  pressures  and  these  also  yield  smaller  discharge 
coefficients. 


54.  The  present  study  on  the  discharge  coefficients  unfortunately  is 
very  incomplete.  It  would  have  been  the  better  procedure  to  base  the  eval¬ 
uation  of  C  on  the  direct  measurement  of  the  mass  rate  of  flow  as  deter¬ 
mined  by  the  simultaneous  temperature  and  pressure  measurements  in  the 
receiving  vessel  as  function  of  time.  It  was  the  intention  to  follow  the 
procedure  in  an  extended  study,  but  due  to  the  lack  of  a  sensitive  tran¬ 
sistor  for  the  temperature  measurements  at  that  time,  the  idea  was  aban¬ 
doned.  Another  precaution  in  the  conduct  of  tests  would  have  been  the 
measurements  of  the  pressures  in  various  locales,  one  important  place  being 
the  vicinity  of  the  orifice.  There  was  only  one  pressure-measuring  element 
used  in  the  tests  made,  and  this  definitely  was  not  sufficient.  Again,  for 
a  more  reliable  control  of  the  tests,  it  would  have  been  desirable  to  con¬ 
nect  the  receiving  vessel  to  a  very  large  discharging  vessel  to  replace  the 
outside  air.  In  this  arrangement  there  would  be  an  additional  means  to 
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determine  the  mass  rate  of  flow  of  the  air. 

55.  The  main  assumption  of  the  above  analysis  leading  to  equation  63 
was  that 

po 

Now  suppose  that  in  an  actual  situation  the  relation  between  the  pressure 
and  the  density  of  the  air  in  the  receiving  vessel  is 

i«  1- 
po  7*o 

where  e  is  a  numerical,  the  value  depending  on  the  circumstances  of  the 
filling;  it  can  be  a  variable,  changing  with  time.  If  the  process  is  to 
approximate  isothermal,  adjustment,  it  may  be  a  number  close  to  7  •  Re¬ 
peating  the  analysis  one  now  obtains,  in  the  place  of  equation  62,  the 
relation 


A/e  &l\ 

dtV 


=  -KCa 


where  K  has  the  same  meaning  as  bef  jre.  When  information  relating  to  e 
is  lacking,  owing  to  the  insufficiency  of  observation,  +he  solution  of 
equation  64  is  uncertain. 

56.  In  any  case  one  may  revert  to  equation  62  or  to  equation  63  as  a 
convention,  and  C  determined  in  this  manner  would  be  referred  to  as  the 
effective  coefficient  of  discharge. 

Pressure  Buildup  in  Pneumatic  Chambers 


57.  In  computing  the  pressure  changes  in  pneumatic  chambers  open  to 
air,  due  consideration  must  be  given  to  the  fact  that  during  the  generation 
of  an  elongated  wave  of  constant  depth  the  air  volume  in  the  chamber  would 
be  changing.  The  change  is  measured  in  terms  of  the  fall  of  the  water  sur¬ 
face  inside.  Ideally,  the  fall  must  be  uniform;  that  is,  if  V  is  the  air 


volume  at  time 


and  Vq  the  initial  value 
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V/y„  =  1  +  n- 

'  u 


(63) 


where  n  is  a  constant.  Simultaneously  the  pressure  buildup  in  the  cham¬ 
ber  must  conform  to  the  expression 


£p/£p0  =  1  -  ^  -  2kgt 


(43  bis) 


especially  in  the  l.ater  stages  of  the  wave  formation.  The  bearing  of  the 
quantities  and  upon  the  wave  height  h  and  the  water  depth  in 
the  channel  has  been  shown  previously  through  equations  51  and  53* 

53.  The  rate  n  here  is  proportional  to  2kg  .  Denoting  the  eleva¬ 
tion  of  the  chamber  ceiling  from  the  channel  bottom  by 


and 


Hence, 


where 


VQ  =  (Hj  -  H  )l 
V  -  (Hj  -  H  +  6)i 


vo  A 


A  *• 


H 


Comparing  this  with  equation.  65 


or 


_  1  d& 
n  A  dt 


n  =  v-j/A 


and  from  equation 


H  -  H, 


or 


n  = 


n  = 


0 


*2k2 
^  -  H  *  2k2 


A 

H-Hq 


(66) 

(67) 


(68) 


(69) 


59 •  In  the  experiments  of  pneumatic  wave  generators  of  the  type 
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considered  here,  the  passage  from  the  pneumatic  chamber  to  the  outside  air 
consisted  either  of  a  short  tube  or  a  rubber  hose  leading  from  the  short 
tube  or  a  circular  orifice.  In  all  these  cases,  once  the  water  was  sucked 
to  level  H  and  the  passage  was  opened  to  air  outside,  the  fall  of  the 
water  surface  inside  and  the  subsequent  increase  of  pressure  in  the  chamber 
took  place  uniformly  without  further  control  of  the  air  passage.  This 
matter  deserves  an  analytical  examination,  and  it  will  be  supposed  that  the 
air  aperture  is  In  the  form  of  an  orifice,  since  the  coefficient  of  dis¬ 
charge  of  this  type  aperture  is  known. 

6a.  We  take  the  general  problem  of  flow  of  air  into  a  vessel  of 
varying  internal  vclume .  The  mass  of  air  in  the  chamber  is  Vp1  and  using 
equation  66 


(70) 


Differentiating  with  respect  to  time 


or 


(71) 


Assuming  an  adiabatic  process  and  remembering  that  the  pressure  inside  the 
chamber  differs  but  little  from  the  outside  atmospheric  pressure. 


Substituting  this  information  in  equation  70, 
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Sinca  Ap  is  small  in  comparison  with  .  the  latter  equation  may  he 
simplified,  ignoring  Ap/Ap^  ,  to 


where 


dm 

dt 


•w.  SH-i) 


d  Ap 
dt  ApQ 


+  KP 


1  d& 
A  dt  _ 


(72) 


N 


P 


(73) 


6l.  It-  will  be  recalled  that  the  mass  rate  of  air  flow  from  outside 
into  a  vessel  will  be  expressed  by 


ff  =  Pq/2  Ca(Ap0)l/2  n/ap/%  (59  bis) 

or  writing 


«  =  Ap/ApQ 


dm  l/2  _  / .  \  l/2  l/2 

dt  =  pO  Ca(V  '  * 


Here,  C  and  a  are  the  coefficient  of  discharge  end  the  area  of  the  ori¬ 
fice,  respectively.  Eliminating  dm/dt  between  equations  59  bis  and  72 
and  making  some  simple  transformations,  the  result  is 


4  ♦!>£♦*,&  I-®** 


(74) 


where 


P0 


V1/2  , 


o(!o_f 

vv 


V'**  “d  K  =  — 


Th's  is  the  general  equation  for  the  buildup  of  tae  press’tre  in  a  vessel 
whose  internal  volume  V  =  V^(l  +  6/A)  changes  with  time.  To  specialize, 
put 

5/A  =  at 

and  introduce  this  in  equation  74;  since  n  will  be  taxen  to  be  a  constant, 
the  result  is 
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.  /lit  1/2 

-(1  t  at]  —  +  N  n  =  KCan  ' 

'  dt-  p 

The  dimensions  of  the  product  KCa  are  the  inverse  of  time. 


/ 

V 


Introduce 


tQ  =  (KCa)'1 
t  =  t/tQ 


With  these,  equation  75  takes  the  dimensionless  form 


-(l  +  mt)  ^  +  mNp  =  s1^ 


(76) 


This  is  the  differential  equation  that  should  apply  for  the  pressure  build¬ 
up  in  the  pneumatic  chamber  when  the  water  surface  is  falling  uniformly. 

The  differential  equation  with  m  a  constant  permits  an  exact  solution. 

For  the  purpose  at  hand,  however,  it  would  be  better  to  consider  a  series 
solution  in.  the  ascending  powers  of  t  . 


it  = 


(77) 


where  the  0  subscript  outside  of  a  parenthesis  indicates  that  the  quan¬ 
tity  is  to  be  evaluated  for  t  =  0  .  Assuming  that  C  is  independent  of 
time,  and  this  would  be  the  case  if  the  coefficient  of  discharge  were  not 
affected  by  the  volume  of  air  in  the  pneumatic  chamber,  one  will  obtain  the 
following  sequence  for  the  coefficients  appearing  in  'auation  77  by  refer¬ 
ring  to  equation  76  and  those  obtained  by  the  successive  differentiation  of 
the  latter  with  time.  The  sequence  is: 


(78) 
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The  juaivtities  in  the  first  parenthetical  term  on  the  right-hand  side  of 
the  last  equation  are  positive,  and  therefore  the  second  time  derivative  of 
«  does  not  vanish.  The  meaning  is  that  the  curve  of  «  versus  time  is 
not  a  straight  line  but  is  curved  slightly.  The  solution  of  equation  76 
does  not  exactly  reproduce  the  linear  variation  of  the  pressure  in  accor¬ 
dance  with  equation  i+3.  It  is,  however,  expected  that  the  time  variation 
of  it  as  given  by  the  second  line  in  equation  78  would  be  sufficiently  cor¬ 
rect.  Identiiying  (dit/dt)^  with  “Skgtg  and  supposing  that  k^  is  small 
in  comparison  with  unity,  the  second  equation  in  equation  78*  after  revert¬ 
ing  to  the  original  variables,  yields 

2k.,  =  N  n  -  KCa 
P 

Introducing  n  from  equation  69,  remembering  that  =  7Pq/Apq  ,  and 
simplifying,  one  has 


KCa 


1  +  7 


H  -  Hq 

Ht-H 


Since  Pq/^Pq  is  a  large  number,  the  above  may  be  written  next  as 


where 


2k2 


KCa 

H  -  Bo  >0 
Sr-Xo  *0 


(79) 


Equati  ">a  79  purports  to  connect  the  rate  of  increase  of  air  pressure  in  the 
pneumatic  chamber  with  the  area  of  the  orifice  opening.  From  the  tests  on 
the  orifices  one  suspects  that  C  is  close  to  unity.  Fairly  reliable 
values  of  C  may  be  obtained,  it  is  believed,  on  the  basis  of  equation  79» 
using  data  of  the  quantities  appearing  in  this  equation.  The  pertinent 
data  are  ;hown  in  table  5,  together  with  the  computed  values  of  C  . 
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For  the  four  tests  the  mean  value  of  C  is  0.95. 

62.  Because  of  the  manner  of  derivation  the  values  of  C  .iusx  re¬ 
ferred  to  apply  for  snail  t  or  small  water  displacement  in  the  pneumatic' 
chamber.  What  the  values  may  be  for  subsequent  times  or  greater  displace¬ 
ment  can  be  obtained  from  equation  74  inserting  in  it  the  observed  values 
of  dn/dt  ,  d&/dt  ,  «  ,  and  6  .  The  values  of  C  that  were  evaluated  in 
this  manner  are  presented  in  table  6.  It  was  mentioned  previously  that 
during  the  action  of  the  pneumatic  generators  the  water  in  the  chamber  re¬ 
mains  immobile  for  a  few  seconds  after  the  orifice  is  opened  to  outside 
air.  Thus,  for  this  initial  period  d&/dt  is  zero.  Accordingly,  a  spe¬ 
cific  value  C  may  be  computed  from  equation  74,  putting  d&/dt  equal  to 
zero.  Values  thus  obtained,  see  figs.  8  and  9>  are  also  included  in 
table  6.  Now  C  is  a  discontinuous  function  of  t  or  6  ,  so  that  if 
equation  74  is  chosen  as  a  basis  for  the  analytical  treatment  of  generator 
action,  this  discontinuity  must  be  borne  in  mind.  This  may  be  accomplished 
in  the  following  manner.  Let  be  the  effective  coefficient  for  the 

initial  peiiod  or  for  the  times  when  d6/dt  vanishes  and  C  the  effective 
coefficient  for  the  later  times.  Denote  the  ratio  of  C  to  by  r 


Insert  this  in  equation  ?4  and  write  in  the  resulting  equation 


Q  =  6/A 

This  gives 

-a ♦  8>  If +  BP I  -  (so) 

Note  that  the  dimensions  of  KG^a  are  that  of  time  inverse;  accordingly, 
put 

tQ  =  (KCa)"1  (81) 


and  introduce  the  dimensionless  time  variable 


T  =  t/t 


0 


in  the  above  equation.  The  result  is 


x  '  dT  p  dt 


(82) 


Now,  when  d^/dr  does  not  vanish,  the  first  term  on  the  left-hand  side  of 
equation  82  is  insignificant  with  respect  to  the  second  term,  is  a 
large  number,  and  tiie  above  equation  for  larger  times  may  be  replaced  by 


(1  ♦  e)  ar  ♦  !!e  !» ,  PA 

'  '  df  r  dr 


(83) 


63-  In  the  example  of  a  numerical  analysis  of  the  pneumatic  gener¬ 
ator  action  the  above  form  of  the  equation  will  be  used  to  represent  the 
flow  of  air  into  the  pneumatic  chamber. 
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PART  V:  INITIAL  CONDITIONS  OF  CHANGING  PRESSURES 

General  Relations  in  Air  and  Water  Flows 


64.  The  discussions  given  in  the  preceding  sections  refer  to  the 
later  stages  of  the  pressure  buildun  in  the  pneumatic  chamber,  after  the 
chamber  has  been  opened  to  the  air  outside.  For  these  later  times  the  in¬ 
crease  of  pressure  is  uniform,  and  the  accompanying  fall  of  water  surface 
in  the  chamber  is  also  uniform.  For  the  initial  times,  just  after  the 
chamber  is  opened  to  air  outside,  the  pressure  change  within  the  chamber 
is  undulatory.  First,  there  is  an  increase  of  pressure  of  considerable 
magnitude,  and  next,  a  pressure  fall.  These  are  followed  by  cyclic  changes 
of  gradually  decreasing  magnitudes. 

65.  Consider  once  more  the  curve  in  fig.  5.  Note  that  the  wave 
does  not  emerge  from  the  chamber  immediately  after  it  is  opened  to  air, 
but  about  0.3  or  0.4  sec  later.  During  this  period  of  time  the  eleva¬ 
tion  of  water  in  the  chamber  remains  the  same;  that  is,  during  this  time 
interval  Vq  is  maintained.  As  a  consequence,  there  is  a  large  increase 
of  pressure  owing  to  entering  air,  over  what  it  would  have  been  if  the 
water  surface  had  started  to  fail  immediately  after  the  chamber  was 
opened.  The  large  increase  of  pressure  forces  the  water  out  of  the 
chamber  almost  suddenly  and  at  a  rate  that  has  the  effect  of  producing 

a  vacuum,  a  depression  of  pressure.  Similar  effects  are  reproduced  in 
the  sequence  but  with  lesser  severity.  The  explanation  is  that  the 
inertia  of  water  is  the  underlying  cause  cf  the  initial  oscillations  of 
the  air  pressure. 

66.  Iu  a  mathematical  treatment  dealing  with  the  initial  condition 
it  i3  necessary  to  determine  the  flow  from  the  chamber  into  the  channel. 

The  expression  for  the  flow  may  be  obtained  from  equation  7a.  In  this, 
write  g(H  -  &)  in  the  place  of  ,  -£p  in  place  of  P  ,  gy  in  place 
of  Gg  ,  and  pg(H^j  +  h  -  y)  in  place  of  P^  .  Assume  that  v^  is 
independent  of  4  and  independent  of  y  .  Thus,  the  energy  equa¬ 
tion  is 
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&  ;y  7- 


Ct'O 


-  ^(Hjj  +  h)  [od(H0  +  h)  +  |  u“]  (Sl») 

I  .'  <ni  •  flow  throug::  tiie  nozzle  is  resisted  because  of  friction,  then  At' 

! to  b<_  replaced  by  At;  +  Ap  where  Ap  ,  a  positive  quantity,  repre- 

lo  1j 

s-'-nts  tne  additional  pressure  difference  to  over  orno  the  resistance.  In- 
t rodu  •  in/  .is  T  and  rememb-  ring-  that  i  =  up(1Iq  +  h)  ,  the  result  may  be 
written 


wh- M  is  a  dimensionless  quantity,  the  value  of  v/hich  depends  on  the 
-■‘ora  of  the  pneumatic  chamber  and  the  nozzle  The  determination  of  it  is 
discussed  in  the  Addendum.  In  view  of  this  latter  relation,  equation  85 
may  no*.:  bo  .-:r itton  as 


=!•«’  “  “  vj  ~  f  «2  +  cs(H  ”  H0)  -  bgh  ~  egb  -  ApL  -  &p 


Dividing  by  eg(H  -  liQ)  or  by  its  equivalent  ApQ  and  remembering  that 
v,  =  do/dt  , 


i‘iJ*  'X  ^ 

c*(!{  -  H  )  2 

'J  ^  o 


2  2 

h  5 

2  g(H  -  HQy  '  H  -  Hq  ‘  H  -  Hq 


1  u?  "  vp 


(37) 


A  s ' I- ■  transformation  will  show  that 


-  -i  )  FT 


"  U0  VW  L  2H0  V-  /  J 


hO 


or  with  a.  sufficient  degree  of  a.pi»‘Ox  '.mat '  on 


2  2  u 

U2  ~  "l  _  1{0  (  h 

iumy "  h  -  v>r 


And  also,  as  it  is  shown  in  the  Addendum, 


ApL 


egorng  -  h  -  n0  ?  VHo 


h£>  i/jl? (b) 

»} 


so  that 


1  U2  ~  V1  ApL  :lO  1 

2  g(H  -  Hq)  og(H  -  H0)  H  -  H0  |_2  2  \  D 


Substituting  these  in  the  energy  equation  the  latter  beconv-o 


fie 


d26 


ew  -  V  dt2 


■  1  - 


H-K0  H  -  Hq  H-HoL 


1  +  i(2° 

2  2  \  D 


)2  _ 


V 


Apo 


(«J) 


Introducing  the  dimensionless  variables 


=  h/H0 
6  =  6/A 


it  =  Ap/&pQ 

t  =  t/tQ 

where  tQ  is  a  characteristic  time,  see  equation  8l,  the  energy  equation 
becomes 


where 


^-§  =  \(l  -  Ag©  -  A.  n  -  Aki)2  -  it)  (89) 


M/A 


;(h  -  H0)t: 
A 


H  -  H. 
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II. 


0 


3  H  -  H. 


and 


.  J5sL.fi +  i(5»f  I 

H  -  Hq  |_2  2  V  D  /  J 


The  relatfon  betwe  n  &  and  h  results  from  the  continuity  condition 


That  i s , 


vl*  =  ’-2(-Ho  +  h) 


S -M>  - 1 -  r(WI:) 


In  terms  of  dimensionless  variables 


d£  , 

dt  "  £A 


(90) 


Inverting 


v-nere 


and 


_  _  .  d6  .  /d Of 
^  A5  di  “  A6  (  dr  ) 


A5  = 


^  •  Vo 


(91) 


^  2 

A6  “  ¥  A5 

67.  Equations  88  and  89  describe  the  flow  cf  xvater  out  of  the  cham¬ 
ber.  To  these  one  must  now  add  the  equation  describing  the  flow  of  air 
Into  the  chamber;  that  is, 


-U  ♦  «>  g  +  *8  g  - 


(92) 


which,  in  conformity  with  equation  83,  yields 


1*2 


and 


I  k 


P  5- 

.r 


I-  r 


i. 


i  ^ 

i  i 

i  [ 

e  t 


'?  f 

f  t 


*e  =  T 


A =  1 


68.  Our  study  of  the  flow  of  air  into  the  pneumatic  chamber  was  not 
sufficiently  complete  to  formulate  a  theoretical  basis  for  determining  Ag 
and  Ay  a  priori.  In  the  numerical  analysis  to  be  shown  subsequently,  the 
values  of  Ag  and  Ay  will  be  inferred  from  te„t  records  t'nat  enable  one 
to  give  n  and  0  as  functions  of  t  . 

69.  Equations  90>  91j  and  92  are  the  basic  relations  describing  *  , 
&  ,  and  q  as  functions  of  t  .  The  desired  solution  may  be  from  the  ap¬ 
propriate  difference  equations  shown  below,  using  the  technique  of 
computers . 

Initial,  conditions: 


Differences:  t  =  nAr 


0  =  o,n*o,it*i,xao 

@L  •<*).•©)“ 

n 

=  g  +(*e) 
n  Vdi/R 

=  q  +  (4^) 

n  \  dT/ 


(93) 


n+l 


At 


‘n+l 


(94) 


Relations  at  t  =  nt; 


&)  -  V1  -  Vn  -  *3\  -  Vg  -  "J 

'  ■  n 

\=A5(i) 

n  n 

(3T)„=-[vi/2-A8(f)n]i<1^-) 


(95) 
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An  Example  of  Numerical  Evaluation 


70.  To  see  whether  the  relations  developed  ir.  the  previous  section 
are  sufficient  to  describe  the  act : on?  of  a  pneumatic  generator,  the  case 
of  Run  g  will  be  considered.  The  graphs  of  the  fall  of  water  surface, 
the  increase  of  internal  pressure,  and  the  height  of  the  issuing  wave  were 
shown  in  fig.  5.  First,  th  characteristic  time  t^  as  defined  by  equa¬ 
tion  Si  needs  to  be  established  by  noting  the  variation  of  pressure  in  the 
•hamber  with  time  prior  to  the  emergence  of  the  wave  from  the  generator, 
i.ie  data  are  shown  in  fig.  12,  and  the  straight  line  drawn  implies  that 

—  KCa  =  1.65  per  sec 

ana  hence 


-  0.3  sec  and  t  =  3.3t 

To  complete  th<-  listing  of  the  pertinent  quantities  of  the  run,  the 
following  ones  are  added: 

H0  =  0.319  ft 

H  =  1.419  ft 

=  1.335  ft 

A  =  0.466  ft 

i  -  2  ft 


The  generator  Is  of  the  elevated  type 
and  M  are  discussed  in  the  Addendum 
peaking  in  equations  89  and  93  are: 

A1  =  2.06 
Ag  =  0.424 
A3  =  0,284 


Pq/pS  =34.1  ft 
Z5p0/pg  =  1.094  ft- 
7  =  1.408 
B  =  66  F 
5  =  1.66 
M  =  6.99 

and  the  appropriate  values  of  t 
With  the  above,  the  constants  ap. 

A^  =  0.682 

a5  =  3.11 

A6  =  7.23 
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71.  Next  must  he  indicated  the  values  of  the  constants  and  Ag 
in  equation  92.  The  connections  to  outside  air  and  to  aspirator  wer<_  made 
through  a  so-called  three-way  valve.  The  characteristics  of  the  air  passage 
are  not  known.  However,  the  limiting  values  of  A.j  and  Ag  can  be  estab¬ 
lished  on  the  basis  of  equation  92  by  inserting  in  it  the  limiting  vs  lues  of 
dfl/dT  and  d 0/ax  .  Hie  data  of  the  observations  are  shown  in  fig.  13  -in 
terms  of  the  dimensionless  variables.  It  is  seen  that  dS/dx  and  dn/dt 
are  practically  constant  for  large  values  of  t  .  For  x  =  10  ,  one  has 

e  =  0.656 
n  =  0.676 
dfi/dx  =  0.0624 
-an/dt  -  0.0254 

Substituting  these  in  equation  94  yields 

0.822  =  O.Ql'21  +  0.0624  Ag 

a  linear  equation  in  the  constants  A^  and  Ag  .  If  it  is  desired  that 
equation  92  yield  the  pressure  increases  observed  during  the  initial  portion 
of  the  limes,  one  must  have 

V1 

and  hence 

Ag=12.7 


72.  It  would  be  instructive  to  inquire  ’f  there  are  differences  be¬ 
tween  the  values  of  Ag  for  the  moderate  x  ,  of  t  and  the  limiting 
value  just  given.  Since  Ag  =  ,  any  variation  in  Ag  would  indicate  a 

corresponding  variation  in  r  :  that  is,  the  ratio  of  the  effective  coef¬ 
ficient  Z  to  its  initial  value  (\  .  Through  equation  90,  expressing 
d^/dx  in  terms  of  tj  ,  one  has  from  equation  92 


*3  = 


=  Ar 


*V2 


+  (l  +  6)drt/dT 


1  +  h) 
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(96) 


The  values  of  Ag  computed  on  this  basis  are  given  In  fig.  1J.  The 
straight  line  drawn  in  the  figure  shows  the  Halting  value  of  Ag  =  12.7. 

The  value  of  Ag  for  moderste_  t  oscillates  about  the  limiting  value. 

This  fact  will  be  ignored,  and  the  numerical  analysis  will  be  attempted 
assuming  tiiat  Ag  is  a  constant  throughout. 

73.  A  Fortran  program  for  the  numerical  solution  of  equations  93, 

94,  and  95  is  shown  in  table  7.  The  numerical  results  are  plotted  in 
fig.  14  and  may  be  compared  with  the  observational  values.  The  initial  and 
final  values  of  ft  are  correctly  computed,  except  that  for  the  moderate 
times  there  is  a  difference  in  the  oscillation  periods.  This  is  due  pri¬ 
marily  to  the  circumstance  that  in  reality  Ag  is  not  a  constant  quantity 
as  assumed  for  the  computation.  The  final  values  of  tj  and  6  are  also 
reproduced  nearly  correctly.  There  ars,  however,  c*arlted  differences  be¬ 
tween  the  observed  and  computed  values- for  the  initial  portion  of  time. 

This  suggests  that  the  relation  between  q  and  dfl/dT  implied  in  equation 
90  or  equation  91  although  valid  for  large  values  of  t  is  not  sufficiently 
valid  for  the  instant  that-  the  wave  is  emerging  from  the  generator  and  ac¬ 
cordingly  for  the  rigor  of  the  analysis  the  necessary  modification  needs  to 
be  introduced. 

?4.  The  above  example  of  the  analysis  shows  that  in  essence  the 
mathematical  formulation  put  forward  is  quite  adequate  to  describe  the 
mechanical  action  of  pneumatic  generators  for  surges,  t.  better  knowledge 
of  the  orifice  or  the  air  passage  characteristics  can  lead  to  a  more  faith¬ 
ful  representation  of  the  actions  transpiring.  The  analysis  employed  by 
2 

Shen,  from  the  point  of  view  of  mechanics,  is  equivalent  to  the  present 
analysis ,  There  sure,  however,  differences  in  the  mathematical  details 
arising  from  the  fact  that  in  establishing  the  flow  equations  Shen  con¬ 
sidered  the  particle  velocity  at  the  nozzle  mouth  whereas  the  present 
analysis  considered  the  wave  height. 
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75.  The  preceding  discussion  of  experiment;;  and  also  of  theory  lias 
shown  that  pneumatic  tanks  would  serve  wcJl  for  the  generation  of  elongated 
surges  moving  with  constant  heights  and  constant  depths.  The  required  con¬ 
trols  and  manipulations  are  truly  simple.  The  channel  waters  a:c  raised  to 
prescribed  heights  in  the  chamber  by  suction,  and  the  chamber  air  suction 
pressure  is  reduced  by  allowing  the  air  from  outside  to  enter  the  chamber 
through  an  aperture,  preferably  a  circular  orifice,  the  aperture  opening 
remaining  the  same  daring  the  filling.  In  this  process  marked  oscillation 
occurs  both  in  the  air  pressure  in  the  chamber  and  in  the  generated  wave, 
due  to  inertia  effect.  This  oscillation  is  of  short  duration,  and  soon 
steady  conditions  are  established.  The  pressure  increases  uniformly  in 
association  with  the  wave  generation  of  the  required  form.  It  is  more  sim¬ 
ple  td  base  the  pneumatic  tank  design  with  reference  to  these  later  condi¬ 
tions.  With  this  understanding  the  relations  needed  for  the  design  are  the 
following: 


h 

«0 

<  Vho 

i)(i  - 

1  h  \  \ 

%  h0;  1  + 

k 

(A) 

2k g 

,  n  4 

’5#  L 
^tv1 

% 

l 

IB) 

^l"1 

+  X  t  \A_  2  ^ 

-O 

1 

(c) 

2L-: 

ml-h 

-  ; 

%  - 
% 

=  *1 

G>) 

**0 

■=  Pg(H  - 

Hq) 

(E) 

KCa 

3=2=  (3 

H  - 

v 

-H  *0) 

(?) 

At 


r»  I  1  /\ ?  wiirf  r.yw  ^  «V»o Kr.Urv*'* 

Zk>  Ruction  pressure  in  chamber  during  the  subsequent  wave  motion, 
poundals/ft^ 

£p.  Suction  pressure  in  chamber  at  the  beginning  of  wave  motion, 

1  poujidals/ft2 

Initial  suction  pressure  in  chamber  to  raise  water  to  height  H, 
poundals/ft^  (Note:  is  a  positive  quantity) 

X  Friction  factor  of  nozzles 

-  o  Density  of  water,  lb/ft^ 

co  Wave  velocity,  ft/sec 

77.  Once  the  design  wave  height  h  is  selected  for  a  given  depth  of 
channel  water  HQ  ,  equation  A  is  the  relation  io  determine  tae  height  H 
to  which  the  water  in  the  chamber  must  be  raised.  For  the  evaluation,  k^ 
is.  assumed,  and  the  resistance  factor  X  is  assigned  a  very  likely  value 
computed  by  the  relation 


discussed  in  the  Addendum.  Here  D  is  the  depth  of  the  narrowest  part  of 
the  deflecting  nozzle.  A  good  selection  for  ^  would  he  a  value  of  0.07 
up  to  .0.10  or  0.12 i 

78.  The  significance  of  k^  is  found  from  equation  D  which  gives 
the  law  of  variation  of  pressure  in  the  pneumatic  chamber  air  space  for  the 
generation  of  an  elongated  wave  of  constant  depth.  The  effective  suetior 
pressure  at  the  start  of  wave  motion  is  £p.  s  a  quantity  less  than  , 

the  suction  pressure  to  raise  the  water  to  a  height  H  .  Equation  E 
gives  the  suction  pressure  in  terms  of  H  and  HQ  .  Nov; 

=  (l  -  k^jAp^  .  The  multiplier  2k^  gives  the  rate  of  cfc  nge  of 
pressure  in  the  pneumatic*  chamber  ail*  space  and  represents  the  quantity 
~  /"^/ApQ  ,  wrere  Zip  is  the  suction  pressure  during  formation  of  the 
v/ave .  Equation  C  gives  the  value  of  the  ratio  2kg/k,  in  terms  of  the 
wave  height  and  the  length  of  the  generator  in  the  longitudinal  direction. 
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79.  Equation  F  shows  that  tire  value  of  2l<^  -s  dejtendent  on  the 
area  of  the  apertures  that  allow  air  to  enter  the  pneumatic  chamber.  Be¬ 
cause  of  this,  relation  one  has  a  certain  latitude  in  selecting  the  value  of 
k1  .  Unfortunately,  relation  F  is  not  insaediately  applicable  to  pneumatic 
wave  generators  of  a  type  where  entrance  of  air  into  the  chamber  is  through 
tortuous  valves  rather  than  through  a  circular  orifice  as  in  some  of  the 
laboratory  tests.  For  such  cases  the  form  of  equation  F,  although  qualita¬ 
tively  correct,  is  not  applicable  quantitatively,  sinee  the  appropriate 
value  of  G  is  not  known.  This  means  that  each  pneumatic  chamber  utiliz¬ 
ing  various  forms  of  valves  needs  to  be  individually  calibrated.  This  will 
cause  no  difficulty  if  in  the  design  of  the  valves  provisions  are  made  to 
change  the  valve  opening  areas  freely  and  arbitrarily. 

BO.  Such  factors  as  particle  velocity,  wave  celerity,  discharges, 
and  discharge  rates  my  be  directly  obtained  from  equations  H,  I,  and  J  if 
desired. 

81.  When  designing  the  pneumatic  generator  to  be  used  in  the  Hilo 
Bay  tsunami  model  the  above  procedure  was  used.  A  sketch  of  the  generator 
as  originally  designed  is  shown  in  fig.  l6.  Assuming  that  the  limiting 
wave  heights  in  the  central  portion  of  the  bay  shall  never  exceed  50  ft, 
and  that  the  depth  of  water  at  the  bay  mouth  is  about  300  ft,  one  has 
h/ELj  =  0.1o7  .  In  the  model  =  1.5  ft,  and  therefore  the  design  wave 
height  is  h  =  0.25  ft  .  The  base  of  the  model  is  80  ft,  and  the  super- 
fieial  area  loOO  ft  .  The  volume  of  the  wave  would  be  5  ftJ,  prorated  per 
foot-  of  length  of  the  mouth. 

82.  The  depth  of  the  nozzle  throat  was  chosen  as  1  ft.  Applying 
equation  ltq.  X  - 0.31  .  lv  was  noted  during  the  experimental  study  that 

:iad  a  value  close  to  0.075  for  the  better  runs  made  with  the  elevated 
generators.  We  shall  accept  this  value.  Placing  =  0.075  ? 
h/IL  =  O.lo?  ,  and  X  =  0.31  ,  equation  A  yields  H  =  6.0  ft  .  This 
is  the  height  the  water  can  be  raised  In  the  chamber.  As  a  reasonable 
value,  !Lg  was  put  equal  to  4.5.  Since  the  maximum  storage  will  be 
(H  -  H^)  and  since  this  equals  5  ft'’,  the  required  chamber  length  'would 
be  1=6  ft  .  The  air  gap  is  arbitrary  and  therefore  was  assigned  the 
value  of  /X  =  2  ft  which  makes  =  8  ft  ,  the  height  of  the  generator.  . 


The  no22le  was  3.5  ft  in  length  and  was  recessed  in  order  to  have  a  compact 
generator.  A  curved  surface  was  introduced  behind  the  nozzle  to  reduce  the 
losses  if  deflection. 

83.  To  know  approximately  what  size  circular  aperture  might  be 

needed  for  the  passage  to  outside  air,  resort  may  be  made  to  equations  F 

and  G.  First,  cne  must  determine  21*w>  .  From  equation  B,  2k^  =  0.071? 

hence  2kg/k^  =  0.95  ,  which  is  in  agreement  with  equation  C.  In  view  of 

the  quantities  Pq/^?0  =  7.5  ,  c Q  =  3120  ft/sec  ,  VQ  =  12  ft^  ,  equation  G 

yields  K  =  363  per  sec  ft**.*  Returning  to  equation  F,  the  required  apor- 

ture  area  is  0.0048  ft  ,  the  value  prorated  per  foot  of  chamber  widths  If 

the  chamber  width  is  10  ft  and  only  one  aperture  is  used,  the  area  of  the 

2 

aperture  would  be  0.045  ft  . 

84.  It  is  discussed  in  the  Addendum  that  the  influence  of  the  air 
jet  from  the  orifice  impinging  upon  the  surface  of  water  would  be  quite 
insignificant,  if  the  air  gap  is  about  ten  times  the  diameter  of  the  cir¬ 
cular  aperture.  For  the  above  value  the  diameter  is  0.24  ft?  thus  the 
air  gap  selected  is  satisfactory. 


—X  2 

*  Another  way  of  expressing  this  woula  be  K  =  0.00276  sec  ft 


ADDENDUM 


This  Addendum  i resent*  a  few  results  of  general  nature,  especially  in 
to  the  treatment  of  flows  (gas  and  liquid),  which  may  prove  to  he 
o-'asslstance  In  clarions  tl.e  numerous  applications  made  in  the  main  text. 

1.  Eulerian  M  of  energy  equation.  The  Eulerian  form  of  energy 
•‘juat'on  (that  is,  the  form  applicable  in  an  area  within  a  fixed  boundary 
„ot  moving  with  the  liquid)  may  be  obtained  from  the  Lagrsngian  form  (that 
is,  the  form  applicable  to  the  same  liquid  portion  in  motion).  Lamb  grves 
the  derivation  for  the  latter.  However,  there  is  some  interest  if  the 
derivation  is  made  in  another  manner.  Consider  the  equations  of  motion  tor 

a  two-dimensional  field: 

du  _  .  &  +  £-i|K  (97) 

at  ■  bx  ciy  bx  Pdx 

dv  .  ^v  +  v®v  =  .^l  -  i  |E  (98) 

dt  +  a^  +  v^  ay  p 

Together  with  these  one  also  has  the  condition  of  continuity 

i  *■£ (0J>  +  i  <cv)  ■ 0  t9?) 


(100) 


Imagine  that  in  the  velocity  field  of  fig.  2b  one  takes  a  closed 
curve  3  fixed  in  space  and  of  area  S  .  Since  the  field  is  two- 
dimensional,  the  flow  is  into  and  out  of  n  cylindrical  surface  which  is 

normal  to  the  plane  of  the  flow  lines.  The  curve  s  is  the  intersection 

of  the  cylinder  with  the  plane  z-0.  Take  the  length  of  the  cylinder  to 
be  unity.  Multiply  equation  97  by  cu  ,  equation  98  by  pv  ,  and  add. 


Hence , 


1  o 

2  0  dt 


q2  tfUP^tivC^  =  -Ouf 


where 


2  2  .  2 
a  -  u  +  v 


(101) 


The  last  equation  may  be  written  also  as 


1  d  2  _  1  c  •  2* 

O  P  £  ~  ~  O  (pUfl  ) 


(svq2)  +4(ss+ari} 
7  2  Vox  oy  / 


-  &  («*>  -  £  <■»“>  *  <ff +  f') 

-a™ -.£(*>  ♦.»(§♦*) 

which  in  view  of  the  continuity  relation,  -equation  100,  simplifies  to 

1  b  2  q2  dp  ^  do  id/  2%  id/  2X 

2  p  <i  +  olr  +  -Q^r  =  “-  'r:  (pu"  )  -  ~  ~7  (pvQ  ) 


2  <5t  eft  2  ox 


^  (puO)  -  ~  (pva) 

i.  (up)  -  i.  (vp)  +  p (—  +  st) 

ox  '  -r/  oy  v  ri  r\ox  oy/ 


ox  oy  '  *  \ox  oy/ 

Multiply  the  two  sides  of  this  equation  by  to  dy  and  integrate  ever  the 
surface  S  delineated  by  the  cnvrve  s  .  £  f  feet  the  integration  using  the 

Green's  theorems 


JUU  +  mV)ds  =  *  J(^  +  S) 


wnere  i  and  m  are  the  directional  cosines  of  the  normal  drawn  inward  to 
the  curve  s  ,  and  U  and  V  each  are  functions  of  x  and  y  .  One  has 

"St  Pq2  +  ®  at  =  f(P  %■  +  "p  +  P  )  (^u  +  mvjds 


J  (l  It  +  r‘  §i )ds  ‘  J~(c  %  *  flc  +  p ) U«  +  rcvjas 


(102) 


If  the  liquid  is  incompressible,  r—  =  0  and  the  above  relation  simplifies 
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fl#  •  ' /< 


*  2  /• 

(  %■  +  Bp)  vnto- +  /  PV 

s  *3 


(103) 


-.-mere 


v  =  iu  +  crv 
n 

This  is  ‘he  Euler i an  form  of  energy  equation  for  incompressible  liquids. 

The  interpretation  is  that  the  rate  of  change  of  kinetic  energy  in  a 
bounded  region  is  equal  to  the  difference  in  the  potential  and  kinetic 
energies  of  the  liquids  entering  and  leaving  the  boundary  and  the  rate  of 
work  done  by  the  pressure  on  the  periphery. 

2.  Internal  energy  of  gases.  Denote  the  internal  energy  >f  gas  per 
unit  mass  by  e  .  Thermodynamics  suggests  two  processes  by  which  the  in¬ 
ternal  energy  of  a  gas  may  be  determined.  In  the  first,  the  gas  is  made  to 
expand  without  the  inflow  or  outflow  of  heat.  Thus, 


Hence, 


p  dv  +  p.e  =  C  ;  v  =  — 


and  utilizing  the  adiabatic  relation 


pp  =  constant 


one  has 


(104) 


a 


6  =  ~  ■—  £  +  e 

7  -  1  P  0 


(105) 


In  the  second  process  heat  is  added  while  the  volume  is  kept  constant. 


Thus, 


dQ  =  -de 


dQ  =  C  dO 
v 


*■'  -*?  „?»  *  -*.■ 


^z-r  •-> 


whore  Cv  is  the  specific  heat  at  constant  volume  and  d&  is  the  change 
in  the  absolute  temperature.  Accordingly 


6  -  V +  e0 


(106) 


For  a  perfect  gas  0  =  p/(Rp)  ,  =  R  +  ,  and  7Cy  =  .  And 

accordingly  the  expressions  in  equations  99  and  100  are  equivalent  to  each 
other.  Directly  from  equation  98 


be  v  bp 
St 

P 

be  =  V  bp 

Sx  “  ^2  ox 
P 

de  _  jg_  dp 
Sy  p2  cty 


(107) 


3.  Energy  equation  for  gases.  Lamb  shows  that  the  energy  equation 


for  gases  is 


ae (I  +  v  + 


w)  =  f 

Js 


p(Ju  +  i»v)ds 


where 


W=/pedS,T=/|  pq2dS  ,  V  =  /  pfl  dS 

•4  *4  *4 

This  is  the  Lagrangian  form.  To  obtain  the  Euleriaa  form  consider  the  last 
integral  in  equation  102.  Using  the  continuity  condition,  equation  100, 

f  i>(^  +  ^)dS  =  ’f  (p^+pU^  +  P  V^)-dS 

•'s  •'s 

and  from  equation  107 

^ -f  pv§)ds 
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f  /  ou  ,  Hv  \  0  _  r  -I  /  -  /  <££  „  £££  \ 

j  p  {  r—  r  t~  I  tj.S  =  -  l  -  to  i  f  1  t  ■*'  r— —  I 

J  \  qy  /  /  o  c  /  \  ox  oy  / 


clS 


- /*  Jr  (p«£)  +  ~  (o/c)  ]  dS 

•4 


Using:  c-he  "•Gnr :  unity  equation,  -«quat  Lon  99,  ar.u  th--  theorem  of  Green 
f  p(|“  +  cy)d£  =  ~  j  (os)dS  +  f  (iu  4-  mv)pe  d* 

*4  *4  *4 

Introducing  this  result  in  sciuatioh  102,  one  has 


i  * tM- 


+  pe  +  o2)dS  =  %*  +  +  €0)  vn^s  +  pv^ds  (l08) 


'S  ~s 

In  this  oft  is  negligible  in  comparison  with  pe  .  Consider  the  expres¬ 
sion 


fi/e 


which  is  the  ratio  of  potential  energy  of  gas  at  a  point  per  unit  mass  to 
the  internal  energy  of  the  same  mass.  As  ft  may  be  measured  from  any  hor- 
isontal  level,  let  the  plane  of  reference  be  the  plane  passing  through  the 
center  of  the  vessel  containing  the  gas  and  let  y^  be  the  maximum  dis¬ 
placement  of  a  gas  molecule  from  this  plane.  Accordingly,  ft  is  gyffl  and 


sy. 


m 


or 


or 


o  sy^p  gy 

§=(,_!) -JL  -  (r  .  D  sbP 


gy. 


7  -  r(7  -  1)  -5s 


'0 


v/heri* 


is  the  velocity  of  sound  for  the  gas.  Since  £>'m/cq  is  a  very 


small  quantity,  it  is  obvious  that  ft  is  negligible  in  comparison  with  £ 
With  this  understanding  equation  108  reduces  to 
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(pq2  +  pe)dS 

which  is  the  E-ulerian  form  of  energy  equation  for  gases.  The  interpreia* 
tion  is  that  the  change  in  the  kinetic  and  internal  energy  of  a  gas  in  a 
fixed  volume  S.l  is  equal  to  the  work  done  by  the  pressures  on  the  boundary 
and  the  difference  in  the  inflow  and  the  outflow  of  internal  and  kinetic 
energies  through  the  boundary  of  area  s.l. 

It  is  shown  in  the  above  that  the  statement  of  energies  is  consistent 
with  the  .equations  of  motion  in  the  absence  of  viscous  dissipation  or  of 
turbulence  and  also  when  there  is  no  heat  flow  through  the  solid  boundaries. 
Similar  derivations  may  be  repeated  for  the  cases  involving  viscous  dissi- 
patioh  dr  turbulence.  However,  for  the  problem  at  hand  it  is  oust  as  well 
that  the  statement  of  the  energies  is  excepted  a  priori. 

4 .  Energy  method  to  determine  the  flow  of  air  into  a  closed  vessel. 
The  previous  derivation  leading  to  equation  6s  was  based  on  the  condition 
of  continuity  of  mass  assuming  a  definite  relation  between  the  density  and 
the  pressure  of  air  in  the  vessel.  -  One  may  examine  the  flow  of  air  into 
the  same  vessel  using  the  method  of  energy.  Biis  may  provide  new  informa¬ 
tion  bn  the  relation  between  density  and  pressure. 

Let  E  be  the  total  energy  of  the  air,  kinetic  and  intrinsic,  con¬ 
tained  in  a  vessel  of  volume  VQ  at  time  t  In  the  absence  of  flow  of 
heat  through  the  solid  boundaries 

It  “  %  +E2  •  (1“) 

where  E^  is  the  rate  at  which  work  is  being  done  or.  the  opening  by  the 
pressures  and  Eg  is  the  rate  of-  convection  of  total  energy  across  the 
boundaries  and  in  this  instance  through  the  opening.  This  is  in  accordance 
with  the  general  result  in  equation  109 .  Let  pQ  and  pQ  be  the  pressure 
and  density  of  air  outside,  p,  and  p1  the  pressure  and  density  inside, 
u^  the  velocity  at  the  entrance,  and  the  cross-sectional  area  of  the 
vena  contracta.  Now, 


=  jT(p  %-  +  pp)  -vrrds  +  J"  pv}>ris  {109) 
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Ej  +  e2  =  u1£-1  /  y-fr  Pr  +  |  p^i  ) 

°r  /  2\ 

r  pi  i  ui  I 
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Assuming  that  the  flow  of  air  into  the  vessel  is  adiabatic 


2  2  7 

°1  =  7  - 


2=i" 

L_fo  x  /1\7 

-if:oL  vpo/  J 


(111) 


. 
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Put  p,  *  -  /ip  and  assume  that  the  square  of"  Zp/pQ  is  negligible  in 

comparison  with  unity,  with  this  approximation  ^  -  :; 

4  -  2&/0Q  1  >■■■-;,-  ;  ' :: - 


a.i  „k&. 
%  y  % 
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Substituting  in  equation  111. 


Bi  *  E=>  =  “iVo  [r  -  i  *  r  -  1  %  *  /(%  )  ] 

Since  (*P/PiJ2  is  negligible  in  comparison  with  unity,  finally 


„  (  7  1  At,  \ 

■h  ^.“Wo(  —  -  —  TJ 


Introducing  the  approximate  value  of  u,  from  equation  112 


Since 


*1  +  S2  ' 


*  Vo  >/f(j4r -^t|) 


C0  '  Wq/Cq 


5i  +  ®2  -•'‘2  °i(  r)  (  r  -  i  '  r  -  1  %)' 


(113) 


The  jet  of  "air  entering  into  the  vess  •  ‘.most  immediately  downstream 
of  the  opening  undergoes  a  turbulent  expansion.  This  circumstance  intro¬ 
duces  an  uncertainty  in  the  estimate  of  internal  energy  inside  if  the  vol¬ 
ume  occupied  by  the  expanding  jet  is  comparable  to  the  volume  of  the  ves¬ 
sel.  On  the  other  hand,  if  the  vessel  volume  is  many  times  larger  than  the 
volume  occupied  by  the  expanding  jet,  the  errors  in  the  energy  estimate 
will  be  small.  The  uncertainty  cores  about  by  one’s  inability  to  account 
properly  for  the  kinetic  energy  in  the  expanding  jet.  This  can  be  broken 
into  two  parts,  the  kinetic  energy  of  the  mean  flow  at  a  point  and  the 
kinetic  energy  of  the  turbulent  fluctuations.  As  the  latter  energy  goes 
into  heat,  one  should  be  concerned  with  the  kinetic  energy  of  the  mean 
flow.  This  will  be  neglected  if  the  volume  V.  of  the  jet  is  a  small 
fraction  of  VQ  ,  the  volume  of  the  vessel .  Thus  E  =  ,  where  E. 

denotes  the  intrinsic  energy  of  the  air  contained  in  the  vessel.  Now 


and  since 


B.  ~  C  6> 
i  v  1 
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and  neglecting  the  small  fraction  in  the  last  set  of  parentheses  of  the 
rig  vt-hana  side 


d 

dt  £pQ 


(116) 


and  this  is  the  same  equation  as  the  one  previously  derived,  equation  62. 
The  implication  is  that  the  densi.ty  relation  assumed  in  the  first  deriva¬ 
tion 
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is  in  accordance  with  the  principle  of  energy  where  the  volume  occupied  by 
the  jet,  V.  ,  is  very  small  in  comparison  with  Vn  ,  the  volume  of  vessel. 

If  the  volume  of  the  receiving  vessel  is  small,  account  must  be  made 
for  the  internal  kinetic  energy.  This  may  be  done  as  follows .  Let  u‘  be 
the  velocity  of  air  at  any  point.  The  internal  kinetic  energy  now  is 


ifv°  2 

sK  =  t/  Pi(u')  dV0 

J  0 


where  dVQ  is  a  volume  element.  If  u^q  is  the  initial  entrance  velocity, 
one  may  express  the  air  velocities  inside  in  terms  of  this  initial  velocity. 
Thus, 


.1  2Vo  T/jUfflo 

-2"loVoj  (uj  VQ 


2 

ignoring  the  small  differences  between  PQ  and  .  Since  u^  =  SAp^/p^, 
a  more  useful  expression  of  is 


ek=0£7~T^o 


(U7) 


vnere 
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Accordingly,  for  the  condition  the  total  internal  energy  is 


vo 

E  *  yrr  (p0  -  +  «®0) 


(U8) 


and,  hence, 


dt  7  -  1  \dt  ^0  dt/ 


(U-9) 
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Proceeding  as  before  and  omiiting  steps  tliat  are  obvious,  one  now  has 


A  &L 

dt  £pQ 


da 

•tt  =  KCa 

US 


(120) 


which  shows  the  effect  of  internal  kinetic  energy  on  the  pressure  changes 
in  the  filling  of  a  receiving  vessel.  The  quantity  a  remaining  unknown 
the  solution  of  the  equation  cannot  be  worked  out.  Speculate  for  a  moment 
that 


a  =  -p£p/ApQ 


(121) 


where  0  is  a  positive  constant.  In  this  it  is  tacitly  taken  that  a  in¬ 
creases  with  time  whereas  Ap  decreases  with  time.  Inserting  this  in  the 
pressure  equation  above,  one  finds 


jL£L-  =  K-£__a*& 

dt  zsp0  *  1  +  P  a  p0 


(122) 


which  accords  with  the  experimental  results  that  the  effective  coefficient 
of  discharge  in  the  cases  of  the  vessel  of  small  volumes  is  less  than  the 
values  obtained  with  larger  volumes. 

A  somewhat  different  situation  arises  when  the  lower  part  of  the 
closed  vessel  contains  water  instead  of  a  rigid  base.  In  the  case  where 
the  surface  of  water  is  close  to  the  air  opening  at  the  top,  the  entering 
air  in  the  form  of  a  jet  would  be  impinging  on  the  water  and  then  would  be 
deflected.  In  this  encounter  work  is  done  against  water  primarily  by  vis¬ 
cous  tractions  and  also  possibly  by  normal  tractions.  Let  the  energy  con¬ 
sumed  in  this  manner  be  of  the  rate  Ev  .  In  our  inability  to  appraise  the 
loss  numerically,  let  it  be  supposed  that  it  is  proportional  to  the  work 
done  by  the  entering  jet  and  the  energy  entering  in;  that  is, 


Ev  =  o2  (£,  +  E2) 

Accordingly,  the  energy  balance  equation  is 

§  -  (i  -  <%Mel  +  V 


(123) 


(124) 
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in  place  of  the  one  shown  by  equation  110.  Repeating  the  analysis  and 
Quitting  the  steps  of  the  transformations,  one  has 


d  Ap  da 

dt  Af>0  dt 


=  KCa(i 


and  adding  that  a  -  ,  finally 


_d  Ap 
dt  ApQ 


(125) 


(126) 


Thus,  in  these  adopted  views  the  effective  coefficient  of  discharge  depends 
on  the  volume  of  the  vessel  and  the  part  of  the  surface  of  the  ’/rater  in  the 
vessel  in  contact  with  the  air  tiet. 

5.  Energy  method  to  determine  the  flow  of  air  into  a  vessel  of 
changing  volumes.  Let  Vq  he  the  volume  at  time  t^  and  V  the  volume 
at  time  t.  V  =  Vq  +  AV  .  Placing  AY  =  Cgb  ,  where  is  the  area  of 

the  lower  moving  surface,  the  velocity  of  the  surface  is  v  =  db/dt  .  In 
agreement  with  equation  109  the  energy  relation  is 


dE 


it  =  <Ei +  EA  -  <Er+  e2>? 


(127) 


where  the  second  term  on  the  right-hard  side  represents  the  rate  of  work 
done  by  the  gas  on  the  moving  surface  and  the  flow  of  energy,  intrinsic  and 

+  Eg  and  ^  have 
the  same  meaning  and  the  same  values  as  in  the  previous  section  except  that 
VQ  will  be  replaced  by  V  . 

The  work  done  by  the  air  in  the  vessel  on  the  moving  surface  is 

*12  P1  dt  2 


kinetic,  associated  with  this  motion.  The  terms  E^ 


and  the  outward  flux  of  energies,  kinetic  and  intrinsic,  is 


Since  the  kinetic  part  is  small  Li  comparison  with  the  intrinsic. 


( 


(e+e)=e  +e  =  — - —  d  ~ 

'  l  2 '2  -,12  +  22  7  -  1  -1  dt 


E  illation  127  now  may  be  written  as 
V  dAp  .  7 


V  dAp  ^  7  dV  „  /Vf  7  77“ 

"  7  -  1  dt  7  -  1  P1  dt  CaC0  \  7  /  7-1  ^ 

Dividing  by  VqAc^  ,  and  writing  V  =  (l  +  ^ ,  one  has 

+  =  KCa 

\  A/  dt  Ap„  A  dt  fL°a  A/Ap^ 


(128) 


N  =  7P0/4P0 
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This  is  the  relation  for  the  rate  of  increase  of  pressure  in  a  vessel  of 
changing  volumes.  It  is  of  the  same  form  as  equation  74,  which  was  derived 
by  assuming  that  the  pressure  and  density  relation  is 

fl  _  ,  1  ^ 

p0  “  X  '  7  p0 

One  may  modify  this  result  according  to  the  conditions  that  the  vol¬ 
ume  of  the  receiving  vessel  is  small  and  the  moving  surface  is  in  the  form 
of  a  water  surface  falling  down.  Introducing  the  corrections  to  take 
account  of  the  internal  kinetic  energy  and  the  rate  of  dissipation 

from  the  surface  contact  Ev  and  repeating  the  steps  of  analysis  as  was 
done  in  the  case  of  a  closed  vessel,  and  omitting  details,  one  obtains 

Assuming  as  before  that  cx  =  -PAp/Ap^ 
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_/1  +  £Yi£E-  +  _i_!E®. 

\  x  A/dt  &pQ  1  +  p  A  dt 


KCa 


In  view  of  the  fact  that  and  P  are  not  known,  it  would  be 
desirable  to  adopt  the  relation 

(128  bis) 

to  describe  the  variation  of  pressure  in  the  receiving  vessel  with  its  vol¬ 
ume  changing  with  time,  with  the  understanding  that  C  ,  the  effective  dis- 
cliarge  coefficient,  could  be  a  function  of  time. 

6.  Turbulent  expansion  of  air  .jets.  It  was  assumed  in  the  develop¬ 
ment  of  the  two  preceding  sections,  that  if  the  volume  of  the  expanding  jet 
is  small  in  comparison  with  the  volume  of  the  vessel  receiving  the  jet, 
then  the  internal  energy  of  the  air  in  the  vessel  is  essentially  intrinsic. 
What  the  jet  size  is,  longitudinally  and  laterally,  will  now  be  discussed. 
It  is  known  that  in  a  one-dimensional  jet  issuing  from  a  slit,  or  in 

an  axial  jet  issuing  from  a  circular  orifice,  the  width  increases  linearly 

9 

with  distance  from  the  opening.  The  result  from  Forthman  on  the  width  2b 
of  the  expanding  jet,  the  jet  issuing  from  a  slit,  is  given  in  fig.  17. 
Axial  distance  x  is  measured  from  a  point  upstream  of  the  entrance  and  at 
a  distance  three  times  the  width  of  the  slit.  Roughly,  now. 


.  (±  +  1)±&L. 
\  A/dt  £{>- 
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b  =  0.158x 


(129) 


The  distributions  of  the  velocities  in  the  various  normal  sections  are  sim¬ 
ilar  to  each  other.  If  u  is  the  longitudinal  component  of  the  velocities 
at  a  point  of  distance  y  from  the  median  plane,  and  u  the  velocity  at 
the  axis  itself,  then 

r  - 

m 

At  the  edges  of  the  jet  u  vanishes  and  also  the  derivative  du/dy  . 
Subject  to  these  conditions,  one  may  adopt  the  expression 
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(130) 
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n  \ 
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as  a  first  approximation  of  the  actual  velocities.  Since  the  variation  of 
b  with  distance  x  is  known,  one  may  utilize  the  velocity  expression  to 


determine  the  value  of  u  on  the  basis  of  momentum.  It  will  be  assumed 

m 


that  the  results  obtained  for  a  one-dimensional  jet  apply  also  to  an  axial 
jet.  Since  the  pressure  is  constant  everywhere,  the  momentum  of  the  liquid 
traversing  a  normal  cross  section  should  equal  the  momentum  of  the 
liquid  entering  the  orifice.  Accordingly, 

b 


h  Jo 


y  cty 


Introducing  the  value  of  the  velocities  from  equation  130 


Hb(!  >-*) 


Since  the  area  of  the  vena  contracta  is  half  the  area  of  the  orifice  and 
u^  is  the  entrance  velocity. 


^0  2  *  2 
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where  d  is  the  diameter  of  the  orifice.  Equating  to  each  other 


and  in  view  of  equation  129 


u 


u. 
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In  the  pneumatic  chamber,  let  the  distance  of  the  water  surface  from 
the  orifice  be  x  .  The  jet  is  deflected  laterally  after  reaching  the 


water  surface.  The  pressure  of  impact  at  the  axial  point  is  denoted  by 


u 


%  =  °i  T 


1  2 
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and,  thus, 


2 

^m  =  1‘13  “  Vi 

x 
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Introducing  from  equation  112,  the  difference  Pq  - 


=  1.13d2/x2 


is  small ,  and 


(131) 


where  Ap  is  the  air  pressure  prevailing  outside  the  jet  and  this  is  also 
tne  pressure  in  the  vena  contracta  for  a  pneumatic  chamber  of  large  dimen¬ 
sions.  The  distance  below  the  orifice  such  that  Aj>a  will  be  equal  to 

O.OlAp  may  be  taken  as  the  effective  length  x.s  of  the  jet.  Thus,  the  . 

1 1 

length  of  the  jet  on  this  basis  is 


And  the  volume  of  the  jet. 


X-  =  10.6d 


Y.  =  31. 2d3 


(132) 


(133) 


If  d  is  1  in. ,  for  example,  the  jet  length  X.  is  10.6  ±.\. ,  and  the  vol- 

*3  J 

ume  of  jet  V.  is  0.018  ft  .  If,  on  the  other  hand,  tip  is  chosen  to 
j  un¬ 

equal  0.0029  Ap  ,  then 


1.  *»  21. 2d 
o 


V.  =  250d"> 
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The  above  results  are  for  an  isothermal ly  flowing  air.  If  ZJp/p0  is 
small  the  results  may  be  applied  also  to  a  case  of  adiabatically  moving  air 
as  in  the  pneumatic  chamber  since  the  errors  in  the  application  can  be 
ignored. 

7.  Kinetic  energy  of  liquid  in  chamber  and  nozzle.  Previously  in 
discussing  the  kinetic  energy  T0,  in  the  pneumatic  chamber  and  the 


nozzle  below,  the  energy  was  expressed  as 


X21  =  I  vi  («  -  Kfi  -  S)i  +  §  u^t  (26  bis) 

in  which  the  first  tern  on  the  right-hand  side  represents  the  kinetic 
•energy  ip  the  upper  part  of  the  chamber,  prismatic:  in  shape,  and  therefore 
of  constant  cross  section,  and  the  second  term  represents  the  kinetic 
energy  in  the  lower  part  of  the  chamber  and  the  nozzle  extension  of  it. 

Let  the  kinetic  energy  of  the  lower  part  be  denoted  by  3^  ;  that  is, 

(m 

where  is  a  numerical  constant,  its  value  being  dependent  on  the  shape 

of  the  nozzle  and  the  type  of  the  pneumatic  generator,  either  low  or  ele¬ 
vated.  For  these  two  types  the  evaluation  of  ii+  is  carried  out 

v 

differently. 

Consider  first  a  low  generator  as  in  fig.  l8a.  The  areas  for  which 
the  kinetic  energy  is  to  be  evaluated  consist  of  the  area  (l)  and  the 

area  (2).  Let  a  denote  the  velocity  at  a  point  in  these  areas, 
q  =  u~  +  v  .  It  suffices  to  write  for  a  point  in  area  (l) 

v  =  v^y/D  and  u  =  v,x/t) 

and  for  a  point  in  area  (2) 

u  =  q^R^/t  and  b  =  D  +  (H^  -  D)x/ft 

where  D  is  the  depth  of  the  nozzle  throat.  the  depth  of  the  nozzle 
mouth,  if  is  the  length  of  the  nozzle  and  the  velocity  at  the  nozzle 
mouth,  How, 

q^dA 

where  dA  is  an  elementary  area  arid  A  is  the  combined  area  of  (l)  and 
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factor  M  decreases  with  &  (that  is,  with  the  fall  of  the  water  surface 
in  the  chamber).  For  the  elevated  chambers  the  term  containing  &  may  be 
ignored.  To  illustrate  this  we  now  evaluate  M  for  the  pneumatic,  chamber 
used  in  the  tests  and  for  the  specific  run  that  was  selected  for 

the  numerical  analysis  previously  discussed.  The  pertinent  data  as  seen  in 
table  1  arei 

HQ  =  0.319  ft  Hg  =  0.872  ft 

H  =  1.419  ft  Hj  =  1.885  ft 

t  ~  2  ft  A  =  0.466  ft 

Equation  139  yields 

M  *  6.99(1  -  0.032  5/A)  (l4o) 

which  shows  that  even  when  the  fall  of  the  water  surface  is  twice  the  ini¬ 
tial  air  gap  the  error  in  neglecting  the  5/A  term  is  only  6  percent. 

8.  Ufozzle  coefficient  of  friction.  In  dealing  with  the  problem  of 
chamber  pressures,  it  became  necessary  to  introduce  in  the  energy  equation 
a  term  to  represent  the  total  loss  of  energy  AL  during  the  time  a  wave 
having  the  length  L  and  the  height  h  is  produced.  As  a  matter  of 
convenience,  the  loss  was  expressed  iii  the  form 

AL  =  Xpgh2^  (52  bis) 

that  is,  as  a  fractional  part  of  the  energy  of  the  wave  produced.  The  pro¬ 
portionality  factor  X  was  referred  to  as  the  nozzle  friction  factor.  If 
the  losses  sire  fron  viscosity,  X  would  be  a  constant  independent  of  the 
wave  height  h  .  However,  as  it  is  readily  understood,  the  losses  origi¬ 
nate  from  turbulence  and  for  this  reason  X  would  be  dependent  also  on 
wave  height.  This  is  a  matter  that  requires  some  elucidation. 

From  above,  since  L  =  cot  ,  the  rate  of  loss  to  be  associated  with 
the  pneumatic  chamber  nozzle  area  is  of  the  form 

H  =  Xpgh^No)  (52a) 


71 


Gin'o  the  loss  in  the  region  arises  from  two  sources,  the  deflection  and 
the  subsequent  expansion,  'cne  rate  of  loss  may  be  expressed  also  as 


AL 

At 


AP.u  D 
h  in 


■where  u  is  velocity  in  the  narrow  part  of  the  throat,  D  the  depth  of 
the  throat,  and  AP^  the  differential  energy  head  associated  with  the  com¬ 
bined  losses.  It  is  customary  to  -write 


AP 

L 


2 

m 


where  r  is  the  loss  coefficient.  In  this  relation  £  is  constant,  its 
value  depending  on  the  form  of  the  passages.  Hence, 


AL 

At 


=  C  §  uh) 


m 


From  the 


condition  of  continuity  u__D  =  , 


(141) 


Comparing  this  with  equation  ?2a, 

2 

^gh2a-=cf(^)  u|H0 


(142) 


Ignoring  secondary  quantities. 


V's  “  to 

and 

2 

u  /gh  •=  h/H& 


and  hence  equation  142  reduces  to 


(143) 
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which  connects  the  loss  coefficient  '  with  the  resistance  factor  X  . 

Table  3  presents  the  values  of  X  for  different  runs  with  the  ele¬ 
vated  tank  together  with  the  wave  height  of  the  run  and  the  water  depth. 

The  individual  values  for  every  run  are  the  averages  from  two  ways  of  eval¬ 
uating  X  .  See  tables  1  and  5.  In  the  last  column  of  table  8  are  given 
the  values  of  £/ 2  reduced  from  X  using  equation  143.  The  average  value 
for  £  is 

5  =  1.66 

which  is  quite  reasonable.  Since  the  nozzle  of  the  generator  used  is  of 
low  flare,  the  expansion  losses  would  be  small  and  the  substantial  part  of 
the  loss  should  be  to  the  flow  deflection  through  the  throat  bend.  In  this 
respect  the  bend  action  is  similar  to  that  in  an  elbow  flow  and  for  which 
the  usually  indicated  coefficient  is  £  =  1.13  • 

An  alternate  expression  for  £PT  in  terms  of  u2  is 


which  may  be  put  in  the  form 


H, 


pg(K  -  Hq)  H  -  Hq 


(144) 


a  result  which  was  used  in  the  section  dealing  with  the  general  relations 
governing  the  flow  of  air  into  the  pneumatic  chamber  and  the  flow  of  water 
out  of  the  chamber. 
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Table  1 

Analysis  of  Elevated  Pneumatic  Generator  Data 

(Chamber  connection  to  air  through  three-way  valve) 
Chajnber  width,  1  ft;  length,  2  ft;  flume  width,  0.995  ft; 

P^Pg  =  31*  ft 


m_  t-  *»  .  ^ 

iauic  a 


V0  = 

Determination 

1  of  Discharge  Coefficient  C 

■3  Pq 

‘‘•M  n  :  5§  - 

34.1  ft,  e  =  66  F,  c0  - 

1125  ft/sec 

C 

v  /v 

^0 

/  0 

d,  in. 

t/a 

sr  - 2-9  n 

-~  =  1  f t 
pg 

Mean 

X  103 

xM 

1.467 

0.33 

0.670 

0.690 

0.680 

109.5 

I.32 

1.002 

0.43 

0.700 

0.711 

0.706 

36.1 

0.90 

0.706 

0.68 

0.722 

o.7ki 

0.732 

34.9 

0.62 

0.608 

0.79 

0.701 

0.739 

0.720 

8.0 

0.54 

0.501 

0.97 

0.714 

0.750 

0.732 

,  4.5 

0.44 

0.437 

1.11 

0.691 

0.711 

0.701 

3.0 

o.4l 

0.360 

1.34 

0.792 

— 

0.792 

1.7 

0.33 

0.257 

1.88 

0.768 

0.833 

0.801 

0.6 

0.23 

Table  3 


Decrease  of  Discharge 

Coefficient  with  Volume 

12- 

pg 

34.1  ft, 

pg 

d  =  1.002  in. 

d  =  0.360  in. 

vo 

Vv0 

VA 

ft3 

c 

Xy'A 

X  103 

c 

x  103 

4.03 

0.720 

0.90 

36.1 

0.820 

0.32 

1.66 

3.05 

0.672 

1.19 

47.8 

0.743 

0.43 

2.20 

2.00 

0.592 

1.81 

73.2 

O.696 

0.66 

3.35 

1.46 

0.551 

2.49 

10.0 

0.564 

0.90 

4.59 

O.92 

0.383 

4.09 

15.9 

0.482 

1.48 

7.00 

Table  4 

Determination  of  Discharge  Coefficient 
from  the  Wave  Data;  Pneumatic  Generator  Tests 

Low  Took  Elevated  Tank 


Run 

A14-7 

Run 

Al4-5 

Run 

Al6-2 

Run 

Ai6-.3._ 

v«> 

0.297 

0.291 

O.298 

0.300 

H(ft) 

1.070 

1.076 

1.471 

1.469 

Vft) 

1.609 

1.609 

1.885 

1.885 

*(ft) 

4.75 

4.75 

2.00 

2.00 

vjft3} 

2.58 

2.53 

0 .825 

0.832 

a^ft2) 

2.73  x  10“3 

1.36  x  10"3 

1.36  X  10‘3 

2.81  X  10“' 

4>0/pg(ft) 

0.771 

0.775 

1.171 

1.168 

K 

3430 

3480 

8670 

8600 

KCa 

3.06 

1.80 

6.1 

13.5 

Ka 

9.28 

4.70 

II.80 

24.20 

C 

0.332 

0.383 

0.516 

0.556 

d(in.) 

0.706 

0.500 

0.500 

0.717 

A(in.) 

2.5 

2.5 

5 

5 

VA 

5.96 

4.22 

2.12 

3.04 

Vv0  x  1q3 

23.4 

7.1 

21.9 

76.8 

Table  S 

Analysis  of  Pneumatic  Generator  Data 
(Chamber  connection  to  air  through  circular  orifice) 
Pq/ps  =  34  ft;  6  =  66  P 


Low  Chamber 

Elevated  Chamber 

Run 

Run 

Run 

Run 

*14-7 

*14-5 

*16-2 

*16-3 

Observed 

h 

0.097 

0.064 

0.066 

0.120 

Ho 

0.297 

0.291 

0.298 

0.300 

h* 

1.070 

1.076 

1.471 

1.469 

1.609 

1.609 

1.885 

1.885 

l 

4.75 

4.75 

2.00 

2.00 

vo 

2.58 

2.58 

0.825 

0.825 

a 

2.73  x  10‘3 

1.36  y  10~3 

1.36  X  10‘3 

2.81  x  10“ 

APf/Pg 

o.m 

0.775 

1.171 

1.168 

^Pj/Pg 

0.680 

0.7^4 

1.081 

0.875 

^>/pg 

0.068 

0.040 

0.112 

0.208 

db/dt 

0.075 

0.045 

0.110 

0.209 

*i 

0.170 

0.090 

0.077 

0.184 

*2 

0.089 

0.052 

0.095 

0.217 

K 

3. *6  x  103 

3.48  x  103  . 

8.67  X  103 

8.60  X  103 

Computed 

X 

0.42 

0.46 

0.46 

0.49 

X 

0.25 

0.08 

0.29 

0.62 

c 

0.85 

1.01 

0.92 

1.01 

Table  6 

Effective  Coefficient  of  Discharge  of  the 


Experimental  Rieumatie  Chambers 


C  - 


Low  Chamber _  Elevated  Chamber 


Run 

Run 

Run 

Run 

t.  sec 

^‘1-7 

CM 

£  . 

*16-3 

0 

0.85 

1.01 

0.92 

1.01 

1 

1.06 

1.19 

1.03 

1.12 

2 

1.14 

1.27 

1.09 

1.39 

3 

1.25 

1.39 

1.17 

1.91 

0* 

0.33 

0.38 

O.52 

0.56 

*  Determined  far  the  initial  instant  of  filling  prior  to  the  issuance  of 
wave  from  chamber.  The  table  indicates  the  discontinuity  in  value  of 
the  effective  discharge  coefficient  in  the  action  of  the  chamber. 
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Fortran  Prognu:1.  foi1  Fneumatit:  Chamber  Analysis 


PRIMP  52 

FORMAT  ( ill  I  ,10X,4;:TIME,  !OX,5HTHFrA,13X,2HPI,12X,3HETA) 

TAU -0.0 

THE-0.0 

DTHE-9.0 

ETA -0.0 

PI---L.0 

READ  50 ,D  El-TAU , FHITAU , A I  ,A2,A3,AM5,A6,A7,A8 

50  FORMAT  (IOF7.3) 

5  FRUIT  51,  TAU ,  THE ,  PI ,  STA 

51  FORMAT  (4FJ5.4) 

IF  (TAU -FTHTAU )  2,4,4 

4  read  5o,deltau,fihtau,ai,ao,a3,a4,a5,a6,a7,a8 

IF  (DSLTAU)  99,99,2 
2  ETA-DTHE*  (A.5-aS*DTHE) 

D2ThE=A  i-^  { 1 . 0-A  ?*THE-A  3*ETA-A4*ETA*ETA-Pl) 

DPI— (A7*SQ‘RTF ( PI ) -A8*DTHE)/{  .  O-FffiE) 

THE=THE-®THE*DELTAU 

PI-PUD  PI*D ELTAU 

DTHE-DTHB+D 2THE*D ELTAU 

TAU-TAU+D ELTAU 

CO  TO  3 

99  STOP 

El© 


TaM»  8 

Resistance  Coefficient  t  for  Elevated  Chamber  Nozzle  Flow 

D  =  C.209  ft 


1  1  •  ■  **»w ,v'TW*Y;  'vwstpiN^i^fr^y'tgj!^ 


Initial  and  subsequent  dispositions 
shallow-water  long  wave 


Elevated  pneumatic  chamber 


ELEVATED  PNEUMATIC 
GENERATOR 


Record  of  chamber  water  elevation,  pressure,  and  wave  height 


I,  SEC 

Fig.  8.  Initial  rise  of  pressure  in  elevated 
pneumatic  chamber;  free  orifice 
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Initial  rise  of  pressure  ir»  a  low  pneumai ic  chamber;  free  orifice 


Discharge  coefficient  as  function  of  vessel 


Fig.  15.  Notation  diagram.  Elevated  pneumatic  generators 


17.  Turbulent  expansion  after  Forthman 


